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CHAPTER 1 

Introduction 



1.1 PROBLEM DEFINITION 



The control of physical systems with a digital computer is becoming more 
and more common. Aircraft autopilots, mass-transit vehicles, oil refineries, 
paper-making machines, and countless electromechanical servomechanisms 
are among the many existing examples. Furthermore, many new digital 
control applications are being stimulated by microprocessor technology in- 
cluding control of various aspects of automobiles and household appliances. 
Among the advantages of digital logic for control are the increased flexibility 
of the control programs and the decision-making or logic capability of digital 
systems, which can be combined with the dynamic control function to meet 
other system requirements, 

The digital controls studied in this book are for closed-loop (feedback) 
systems in which the dynamic response of the process being controlled is a 
major consideration in the design. A typical topology of the elementary type 
of system that will occupy most of our attention is sketched schematically in 
Fig. 1.1. This figure will help to define our basic notation and to hitroduce 
^se^^j&iaii fejEfa^ ^ 

^^ ^^^ffi ilS^i Pf ^ ^; — ^ ^^ ^^^ ^ controllM B cMleo^re^ 

response requires control action. 

By "satisfactory response'' we mean that the plant output, y{t), is to be 
forced to follow or trade the reference input, r(f), despite the presence of 
disturbance inputs to the plant [w{t) in Fig. 1.1] and despite errors m the 
sensor [represented by v{t) in Fig. 1.1]. It is also essential that the tracking 
succeed even if the dynamics of the plant shoidd change somewhat during 
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Notation: 

r - reference or command inputs 
u ~ control or actuator input signal 
y = controlled or output signal 

y = instrument or sensor output, usually an approximation to or estimate 
ofy, (For any variable, say 0, the notation 8 is now commonly taken 
from statistics to mean an estimate of $.) 

i = r-y = indicated error 

e = r-y = system error 

w = disturbance input to the plant 

V = disturbance or noise in the sensor 

A/D = analog-to^igital converter 

D/A = digital-to-analog converter 

Figure 1.1 Block diagram of a basic control system. 

the operation. The process of holding y{t) close to r(<), including the case 
where r = 0, is referred to generally as the process of regvlation. A system 
that has good regulation in the presence of disturbance signals is said to 
have good disturbance rejection, A system that has good regulation in the 
face of changes in the plant parameters is said to have low sensitivity to 
these parameters. A system that has both good disturbance rejection and 
low sensitivity we call robust. 

The means by which robust regulation is to be accompUshed is through 
the control inputs to fehe plant [ti(t) in Fig. 1.1], It was discovered long ago^ 
that a scheme of feedback wherem the plant output is measured (or sensed) 
and compared directly with the reference input has many advantages m the 
effort to design robust controls over systems that do not use such feedback. 
Much of our effort in later parts of this book wiU be devoted to iUustrating 
this discovery and demonstratmg how to exploit the advantages of feedback. 
However, the problem of control as discussed thus far is m no way restricted 



^See especially the book by Bode (1945). 
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Figure 1.2 Plot of output versus input characteristics of the A/D converter. 

a digital signal Not surprisingly, digital computers in this book process 
digital signals. 

In a real sense the problems of analysis and design of digital controls 
are concerned with taking account of the effects of the sampling period T 
and the quantization size q. If both T and q are extremely small (sampling 
frequency 50 or more times the system bandwidth with a 16-bit word size), 
digital signals are nearly continuous, and continuous methods of analysis and 
design can be used. The resulting design could be converted to the digital 
format for implementation by using the emulation method described in 
Chapter 5. We will be interested in this text in gaining an understanding of 
the effects of all sample rates, fast and slow, and the effects of quantization 
for large and small word sizes. Many systems are originally conceived with 
fast sample rates, and the computer is specified and frozen early in the 
design cycle; however, as the designs evolve, more demands are placed on the 
system, and the only way to acconunodate the increased computer load is to 
slow down the sample rate. Furthermore, for cost-sensitive digital systems, 
thej|b£&t|de^ the 
required jo^TTH^t'^tT^ 
anl^tfersm^U^ttwordisizei^We^^^ 

T and q separately. We first consider q to be zero and study discrete and 
sampled-data (combined discrete and continuous) systems that are linear. 
In Chapter 7 we will analyze in more detail the source and the effects of 
quantization, and we will discuss in Chapters 5 and 10 specific effects of 
sample-rate selection. 

It is worthy to note that the single most important impact of imple- 
menting a control system digitally is the delay associated with the D/A 
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Figure 1.3 The delay due to the hold operation. 
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converter. Each value of u{kT) in Fig. 1.1 is typically held constant^ until 
the next value is available from the computer. Thus the continuous value of 
u{t) consists of steps (see Fig. 1.3) that, on the average, lag u{kT) by r/2, 
as shown by the dashed line in the figure. If one simply incorporates ttds T/2 
delay in a continuous analysis of a digital system, excellent agreement results 
for many reasonable sample rates. This point will be explained further in 
Chapters 3 and 5. 

Our approach to the design of digital controls is to assume a background 
in continuous systen[is and to relate the comparable digital problem to its 
continuous counterpart. We will develop the essential results, from the be- 
ginning, in the domain of discrete systems, but we will call upon previous 
experience in continuous-system analysis and in design to give alternative 
viewpoints and deeper understanding of the results. In order to make mean- 
ingful these references to a background in continuous-system design, we will 
review the concepts and define our notation as required. 

1.2 EXAMPLE SYSTEMS FOR STUDY 

In order to guide the discussion in the following chapters we have developed 
models for six example control problems in Appendix A. The first of these 
is a satellite attitude control problem in which the plant transfer function is 
the double integrator 



Gi{s) = 



(1.1) 



This example is simple, but with two poles on the stability boundary, this 
transfer frmction must be controlled with care. The second example is a 



^CaUed a Zero Order Hold or ZOH. 
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servomechanism, motivated by a radar tracking antenna- The basic control 
dynamics of the amplifier-motor-load system is second order and has a 
transfer function with one pole at the origin and one pole on the negative 
real axis. The transfer function can be normalized to be 



G2{s) = l/s[is + l)l 



(1.2) 



The third example to be used to illustrate digital control comes from the 
process industries — ^a mixing process. The normalized transfer function that 
represents this system is 



(1.3) 



The obvious feature of this model is the delay term, e"^-**, which is called 
a transportation delay. Such delays are very difficult to control, and the 
example affords us an opportimity to study these difficulties in the con- 
text of digital control. Also, because there is no mtemal integrator in this 
model, good steady-state accuracy requires integral control, which presents 
an interesting but conunon design problem. 

In many electromechanical systems, there is a flexible mechanical mem- 
ber between the actuator and the position sensor. The effect of this flexibility 
is to introduce a complex pole into the transfer function that often has very 
low damping and creates difficulties in providing a fast, stable system. This 
situation arises frequently enough to warrant a name: the noncolocated sen- 
sor and actuator. An example is the control of a robot arm where the end 
pomt of the arm is measured directly. The normalized transfer function for 
this kind of system is 



G4(S) = 



(1.4) 



A variation of this problem occurs if the sensor is on the same mass as the 
actuator, the colocated case. This situation can occur in the positioning of a 
large radio telescope or in the control of a satellite with flexible solar panels. 
In this case the transfer function includes complex zeros with ahnost no 
damping in addition to the Ughtly damped complex poles* The normalized 
transfer function is an extension of (1.4) and is given by 



Gs{s) = 



(1.5) 



Our sixth example il 
it has two controls and 
comes from the process 
and total pressure (hea 
methods are the most ( 
also be described by tb 



where A{s) is the systc 
A seventh example 
magnetic tape drive wl 
of mterest are the tap 
(1.6). 



1.3 OVERVIEW « 

An overview of the pat 
will be useful before 
we place systems of in 
the signals present. Tl 
digital systems. 

In discrete systems 
with an analysis of thej 
signals and "pulse" -tr 
We also develop discr 
sampled, systems tha 
equations and give exs 
descriptions. Having t 
the dynamic response 

A samf^bsd?data s 
often it is important 
For example, with a i 
tween sample instants 
concerned with the qv 
as they might emerge 
necessary for providii 
action typically occui 
to data extrapolation 



1.3 OVERVIEW OF DESIGN APPROACH 



e basic control 
der and has a 
n the negative 



(1.2) 

omes from the 
r function that 



(1.3) 

v^hich is called 
ntrol, and the 
es in the con- 
egrator in this 
v^hich presents 

^ chanical mem- 
'this flexibility 

. often has very 
e system. Thfe 
ncolocated sen- 
where the end 
er function for 



(1.4) 

le solar panels, 
ith almost no 
he normalized 



(1.5) 



Our sixth example illustrates mtiltivariable control. Its key feature is that 
it has two controls and two outputs, so the transfer function is a matrix. It 
comes from the process industry and illustrates control of the liquid level 
and total pressure (head) in the head box of a paper machine. State space 
methods are the most convenient to describe this example; however, it can 
also be described by the transfer function 



A(a)lc(a) d{$) 



6(.)]^ 



(1.6) 



where A(^) is the system characteristic equation. 

A seventh example, used in a multivariable design in Chapter 9, is a 
magnetic tape drive where the tape is driven at both ends and the outputs 
of interest are the tape position and tension. It also can be described by 
(L6). 



1.3 OVERVIEW OF DESIGN APPROACH 

An overview of the path we plan to take toward the design of digital controls 
will be useful before we begin the specific details. As mentioned above, 
we place systems of interest in three categories according to the nature of 
the signals present. These are discrete systems, sampled-data systems, and 
digital sjrstems* 

In discrete sj^tems all signals vary at discrete times only. We will begin 
with an analysis of these in Chapter 2 and develop the transform of discrete 
signals and "pulses-transfer functions for linear constant discrete systems. 
We also develop discrete transfer functions of continuous systems that are 
sampled, systems that are called sampled-data systems. We develop the 
equations and give examples using both transform methods and state-space 
descriptions. Having the discrete transfer functions, we consider the issue of 
the dyna mic respon se of discrete systems, 

often it js ^portant to be able to compute the continuous time r^poiiie. 
Fbr exampfe,^^ii^ lo^ s^^ fi&ere cM be^ipifi^a^^^ 
tween sample instants. Such situations are studied in Chapter 3. Here we are 
concerned with the question of data extrapolation to convert discrete signals 
as they might emerge from a digital computer into the continuous signals 
necessary for providing the input to one of the plants described above. This 
action tsrpically occurs in conjunction with the D/A conversion. In addition 
to data extrapolation, we consider the analysis of sampled signals from the 
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viewpoint of continuous analysis. For this purpose we introduce impulse 
modulation as a model of sampling, and we use Fourier analysis to give 
a clear picture for the ambiguity that can arise between continuous and 
discrete signals, also known as aliasing. The plain fact is that more than 
one continuous signal can result m exactly the same sample values. If a 
sinusoidal signal, yi, at frequency fi has the same samples as a sinusoid j^j 
of a different frequency /a, yi is said to be an alias of j/2- A corollary of 
aUasing is the sampling theorem, which specifies the conditions necessary if 
this ambiguity is to be removed and only one contmuous signal aUowed to 
correspond to a given set of samples. 

As a special case of discrete systems and as the basis for the emulation de- 
sign method, we consider discrete equivalents to continuous systems, which 
w one aspect of the field of digital filters. Digital filters are discrete systems 
designed to process discrete signals in such a fashion that the digital device 
(a digital computer, for example) can be used to replace a continuous filter. 
Our treatment in Chapter 4 will concentrate on the use of discrete filtering 
techniques to find discrete equivalents of continuous-control compensator 
transfer functions. Again, both transform methods and state-space methods 
are developed to help understanding and computation of particular cases of 
interest. 

Once we have developed the tools of analysis for discrete and sampled 
systems we can begin the design of feedback controls. Here we divide our 
techniques into two categories: transform^ and state-space* methods. In 
Chapter 5 we study the transform methods of the root locus and the fre- 
quenqr response as they can be used to design digital control systems. The 
use of state-space techniques for design is introduced in Chapter 6 For 
purposes of understandmg the design method, we rely mainly on "pole as- 
signment," a scheme for forcmg the closed-loop poles to be in desirable 
locations. We discuss the selection of the desired pole locations and point 
out the advantages of using the optimal control methods covered in Chap- 
ter 9. Chapter 6 includes control design using feedback of aU the "state 
variables" as well as methods for estimating the state variables that do not 
have sensors directly on them. A study of quantization effects m Chapter 
7 both presents a "worse-case" analysis and introduces the idea of random 

^Named because they use the Laplace or Fourier transform to represent systems. 

lUe state space is an extension of the space of displacement and velocity used in 
physics. Mudi that is called modem control theory uses diflFerential equations in 
state-space form. We mtroduce this representation in Chapter 2 and use it exten- 
sively afterwards, especially in Chapters 6 and 9. 
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1.4 COMPUTER-AroED DESIGN 9 



signals in order to describe a method for treating the "average" effects of 
this important nonlinearity. 

The last five chapters cover more advanced topics that are essential for 
most complete designs. The first of these topics is identification, introduced 
in Chapter 8. Here the matter of model making is extended to the use of 
experimental data to verify and correct a theoretical model or to supply a 
dynamic description based only on input-output data. Only the most ele- 
mentary of the concepts in this enormous field can be covered, of course. We 
present the method of least squares and some of the concepts of maximum 
likelihood. In Chapter 9 the topic of optunal control is introduced, with em- 
phasis on the steady-state solution for linear constant discrete systems with 
quadratic loss functions. The results are a valuable part of the designer's 
repertoire and are the only techniques presented here suitable for handling 
multivariable designs. Chapter 10 presents a topic with specific application 
to sampled-data and digital controls: the question of sampUng-rate selec- 
tion. In our earlier analysis we develop methods for examining the effects of 
different sampling rates, but in this chapter we consider for the first time 
the question of sample rate as a design parameter. 

In Chapter 11, an introduction to the most important issues and tech- 
niques for the analysis and design of nonlinear sampled-data systems is 
given. The analysis methods treated are the describing function, equivalent 
Unearization, and Lyapunov's second method of stability analysis. Design 
techniques described are the use of mverse nonlinearity, optimal control (es- 
pecially time-optimal control), and adaptive control. In Chapter 12, many 
generic issues concerned with the practical appUcation of digital controls 
are described. The chapter includes a case study of a disk-drive design, and 
treatment of both implementation and manufacturing issues are discussed. 



1.4 COMPUTER-AIDED DESIGN 

As with any engineering design method, design of control systems requires 
many computations that are greatly Jacifitated by a good Ubrary of well- 
documented computer programs. In designing practical digital control sys- 
tems, and especially in iterating through the methods many times to meet 
essential specifications, an interactive computer-juded design (CAD) package 
with simple access to plotting graphics is crucial. Many commercial control 
system CAD packages are available which satisfy that need, and much of the 
discussion in the book assumes that a designer has access to one of them. 
Specific CAD routines that can be used for performing calculations for a 
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topic in the book are indicated in the text, and the routine names for some 
of the popular CAD packages are listed in Table E.l in Appendix E 

CAD support for a designer is universal; however, it is essential that the 
designer is able to work out very simple problems by hand in order to have 
some idea about the reasonableness of the computer's answers 



1.5 SUMMARY 

In this chapter we introduced the subject of digital control design and de- 
fined the variables that characterize discrete, sampled-data, and digital con- 
trol systems. Pointing out that we will be mamly concerned with sunple 
smgle-loop hnear constant systems, we gave the (highly simplified) transfer 
functions of six physical examples, which will be used to illustrate the analy- 
sis and design techniques as they are presented m the book. These examples 
represent aspects of a space vehicle, a servomechanism, temperature control 
a flexible structure, a paper machine and a magnetic tape drive 

Fmally, we gave an overview of our study, from the analysis of discrete 
and sampled-data systems to the design of digital controls via transform 
and state-space techniques. We pomted out that, foUowing development of 
design methods for the most elementary case, we will consider the topics 
of quantization, identification, Imear optunal control, selection of sampling 
rate, nonhnear controls, and issues m the unplementation of digital control 

1.6 SUGGESTIONS FOR FURTHER 
READING 

Several histori^ of feedback control are readUy available, mcluding a Scien- 
tific Amencan Book (1955), and the study of Mayr (1970). A good ^cussion 
Of the histonca^ developments of control is given by Dorf (1980) and by 
Fortmann and Hitz (1977), and many other references are cited by theL 
authore for the mterested reader. One of the earUest published studies of 
control systems operating on discrete time data (sampled-data systems in 
our termmoloar) IS given by Hurewicz in Chapter 5 of the book by James, 
Nichols, and Phillips (1947). ' 

nh J?* '^^f of trackmg and robustness embody many elements of the 
^ectives of control system design. The concept of tracking contains the 
r^ments of system stability, good transient response, and good steady- 
state accuracy, aU concepts fundamental to every control system. Robustn«s 
IS a property essential to good performance in practical designs because real 
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narameters are subject to change and because external, unwanted signals 
invade every system. Discussion of performance specification of control s^^^ 
ems is given in most books on introductory control, mcluding Frankhn, 
pTell, Sid Emami-Naeini (1986). We will study these matters m later 
chapters with particular reference to digital control design 

A comprehensive treatment of conversion techmques between analog and 
dirital signals is given by Hnatek (1976). . r 

A comprehensive text concerned with writing equations of motion for 
systems in a form suitable for control studies is Cannon (1967). 



PROBLEMS AND EXERCISES 

1.1 Suppose a radar search antemia at San FVancfeco fP^^J'^'^'^ ^^^^^ 
and data points corresponding to the position of flight 1081 are plotted on the 
SlXrscreen once^r antemia revolution. Fhght 1081 is travehng chrectly 
Sward the airport at 540 mi/hr. A feedback control s^^tem is established tljough 
Srontroller who gives course corrections to the pilot He wishes to do so each 9^ 
of travel of the aircraft, and his instructions consist of course headmgs m mtegral 

^'^'^^^t is the sampUng rate, in seconds, of the range signal plotted on the 

radar screen? , . ^ i.- o 

b) What is the sampling rate in seconds, of the controUer's instructions? 

c) Identify the foUowing signals as continuous, discrete, or digital: 

i) the aircraft's range from the airport, 

ii) the range data as plotted on the radar screen, 
ni) the controUer's instructions to the pilot, 

iv) the pact's actions on the aircraft control surfaces. 

d) Is this a continuous, sampled-data, or digital control system? 

e Show that it is possible for the pilot of flight 1081 to fly a -^^^ 
which would show up as a straight line on the controller's screen^ What ,s 
the "(t^t) W^^ey of a Sifl'isoidal zigzag, course which wiU be hidden 

from the controller's radar? 
1^ a) Design a continuom lead compensation for the satemte attitude control 
example described by (1.1) so that the complex roots are at s = -4.4±j4.4 

rad/ sec 

b) Assuming the compensation is to be implemented digitally, approximate 
the lag from the ZOH to be 

2/T 



ftl 



Ghis) = Jljlj/T' 
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and determine the root locations for sample rates, a;« = 5 Hz, 10 Hz, and 
20 Hz where T = l/a;, sec if u;, is in units of Hz. 
c) The closed-loop system with roots as specified by (a) has a rise time of 
about 0.3 sec. How fast do you think one should sample in order to have 
a reasonably smooth response? 

1,3 Convert a continuous compensation, 

to discrete hy writing D(s) as a differential equation and then using the definition 
of a derivative, 



x= 11m 
r-»o 



x(k)-x{k-l) 



CHAPTI 

Linec 
Anab 



to convert the differential equation to an equation that determines u{k) given 
u{k - 1), e{k), and e{k - 1). Assume that T is small enough so that the equation 
for X is approximately correct. 

1.4 If a signal varies between 0 and 10 volts (called the dynamic range) and it is 
required that the signal must be represented in the digital computer to the nearest 
5 millivolts, that is, if the resolution must be 5 mv, determine how many bits the 
analog-to-digital convertor must have. 
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Linear, Discrete, Dynamic-Systems 
Analysis: The ^TVansform 
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2.1 INTRODUCTION 

The unique element in the structiure of Fig. 1.1 is the digital computer. The 
fundamental charpcter of the digital computer is that it takes a finite time 
to compute answers, and it does so with only finite precision. The purpose 
of this chapter is to develop tools of analysis necessary to imderstand and 
to guide the design of programs for a computer acting as a linear, dynamic 
control component. Needless to say, digital computers can do many thinigs 
other than control linear, dynamic systems; it is our purpose in this chapter 
to examine their char£u:teristics when doing this elementary control task and 
to develop the basic analysis tools needed to write programs for a real-time 
control computer. 



2>2 Lja>fEAR D IFFERENCE EQUATIONS 

gl^6fatlfe¥8iimaawM ^^ to the coinRliter 

80 that y{kT) = y{kT). The job of the cor|^B^ 
pie values and compute in some fE^hion the signals to be put out through 
the digital-to-analog converter (D/A). The characteristics of the A/D and 
D/A converters will be discussed later. Here we consider the treatment 
of the data inside the computer. Suppose we call the input signals up to 
the kih sample eo,ei,e2, . . . , ejt, and the output signals prior to that time 
t^it ti2f • • • 9 tijb.i. Then, to get the next output, we have the machine com- 
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pute some function, which we can express in symbolic fonn as 

«/fc = /(eo, • • • , ejt; uo, . . . , «jt-i). 



(2.1) 



Because we plan to emphasize the elementary and the dynamic possibilities, 
we assume that the function / in (2.1) is linear and depends on only a finite 
nimiber of past e's and u's. Thus we write 



Uk = -OlUjfe_i-02«fc-2 (^n^k-n+boek + bxek^i + '^' + bmek-m^ (2.2) 

Equation (2.2) is called a linear recurrence equation or difference equation 
and, as we shall see, has many similarities with a linear differential equation. 
The name "difference equation" derives from the fact that we could write 
(2.2) using Uk plus the differences in u^j which are defined as 



(first difference), 
(second difference), 
(nth difference). 



(2.3) 



If we solve (2.3) for the values of Uk,Uk-u and in terms of differ- 
ences, we find 

«it = tiik, 

Thus, for a second-order equation with coeflScients oi,a2, and bo (we let 
6i = &2 = 0 for simphcity), we find the equivalent difference equation to be 

a2V^Uk - (ai + 2a2)Vuk + (02 + ai + l)uk = feoe*. 

Although the two forms are equivalent, the recurrence form of (2.2) is more 
com^nieirth ior^ <50^^ mSl drop ^^^^^g^Q^p^i^i^ 

ferences. We will continue, however, to refer to our equations as "differ- 
ence equations." If the a's and b's in (2,2) are constant, «hen tli% Gdmjjuter 
is solving a constant-coefficient difference equation (CODE). We plan to 
demonstrate later that with such equations the computer can control linear 
constant dynamic systems and approximate most of the other tasks of Imear, 
constant, dynamic systems, including performing the functions of electronic 
filters. To do so, it is necessary first to examine methods of obtaining solu- 
tions to (2.2) and to study the general properties of these solutions. 
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Figure 2.1 The Fibonacci numbers. 

To solve a specific CODE is an elementary matter. We need a starting 
time (A;- value) and some initial conditions to characterize the contents of the 
computer memory at this time. For example, suppose we take the case 



(2.4) 



and start at A; = 2. Here there are no input values, and to compute U2 we 
need to know the (initial) values for tio and ui. Let us take them to be 
tiQ = ui = 1. The first nine values are 1,1,2,3,5,8,13,21,34 — A plot of 
the values of Uk versus k is shown in Fig. 2.1. 

The results, the Fibonacci numbers, are named after the thirteenth- 
century mathematician^ who studied them. For example, (2.4) has been 
used to model the growth of rabbits in a protected environment.^ However 
that may be, the output of the system represented by (2.4) would seem to be 



^Leonardo Fibonacci of Pisa, who introduced Arabic notation to the Latin world 
about 1200 A.D. 

^ Wilde (1964). Assume that Uk represents pairs of rabbits and that babies are bom 
in pairs. Assume that no rabbits die and that a new pair begin reproduction after 
one period. Thus at time we have all the old rabbits, ujb-i, plus the newborn 
pairs born to the matiue rabbits, which are 
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growing, to say the least. If the response of a dynamic system to any finite 
initial conditions can grow without bomid, we call the system unstable. We 
would like to be able to examine equations like (2,2) and, without having to 
solve them explicitly, see if they are stable or unstable and even understand 
the general shape of the solution. 

One approach to solvmg this problem is to assume a form for the solution 
with unknown constants and to solve for the constants to match the given 
initial conditions. For continuous, ordinary, diiflferential equations that are 
constant and linear, exponential solutions of the form are used. In the 
case of linear, constant, difference equations, it turns out that solutions of 
the form will do where z has the role of s and A; is the discrete independent 
variable replacmg time, t Consider (2.4). If we assume that u{k) = Az^, we 
get the equation 

Now if we assume z ^ 0 and A 0, we can divide by A and multiply by 
z *, with the result 

1 = ^-1 + ^-2 



or 



z^^z + 1. 



This polynomial of second degree has two solutions, z = 1/2 ± ^5/2. Let's 
call these zi and Z2. Since our equation is linear, a sum of the individual 
solutions will also be a solution. Thus, we have found that a solution to (2.4) 
is of the form 

tl^) = AiZi''+A2Z2^. 

We can solve for the unknown constants by requiring that this general solu- 
tion satisfy the specific initial conditions given. If we substitute ib = 0 and 
= 1, we obtain the simultaneous equations 
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These equations are easily solved^ to give 

1 + V5 



A2 = 



2>/5 ' 

VE-1 

2>/5 



And now we have the complete solution of (2.4) in a closed form. Further- 
more, we can see that since = (1 + \/5)/2 is greater than 1, the term in 
zi^ will grow without boimd as k grows, which confinhs our suspicion that 
the equation represents an imstable system. We can generalize this result. 
The equation in z that we obtain after we substitute u = is a polynomial 
in z known as the characteristic equation of the difference equation. If any 
solution of this equation is outside the unit curcle (has a magnitude greater 
than one), this corresponding difference equation is unstable in the specific 
sense that for some finite initial conditions the solution will grow without 
boimd as time goes to infinity. If all the roots of the characteristic equation 
are inside the unit circle, the corresponding difference equation is stable. 

Example 2.1: Is the equation 

u{k) = 0.9u{k - 1) ~ 0.2w(fc - 2) 

stable? The characteristic equation is 

z^ - 0.9z + 0.2 = 0, 

and the characteristic roots are 2 = 0.5 and z = 0.4. Since both these 
roots are inside the unit circle, the equation is stable. 

As an ex2uaq]|g bf lW o^iM'W a diffi^M^ aiT^ieMail 
input, we consider the discrete approximation to integration. Suppose we 
have a continuous signal, e(t), of which a segment is sketched in Fig. 2.2, 
and we wish to compute an approximation to the integral 



= f'eit)dt, 
Jo 



(2.5) 



'Do it! 
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(2.6) 



we are led to a simple formula for discrete (trapezoid rule) integration: 

T 

»k = Uk-i + -{ek + ek-i). (2.7) 
If e(t) = then Ck = kT and substitution of Uk = (T^/2)k^ satisfJes (9 7\ 
ht i '^T** f '^l'^'*'^' J "'^^ approximate the area underThe cu^e 

A third possibiUty is the Backward Rwtangular Rule, given by 
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2.3 THE DISCRETE TRANSFER FUNCTION 



2.3.1 The z-TVansform 

If a signal has discrete values cq, ci , . 
the signal as the function*'* 



I e*, . . . we define the z-transfonn of 



-ifc 



(2.8) 



signals thTare zernor 'TT^^^^^ ^^e^z-". For 

take the one-sided tr^fo™ of t T '^^""'^^^ To 

problems. We select the tlStfrT* ^ t"*^^ transform are invited by the 
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and we assume we can find values of ro and i?b as bounds on the magnitude 
of the complex variable z for which the series (2.8) converges. A discussion 
of convergence is deferred until Section 2.7. 



Example 2.2: As an example to iUustrate (2.8), consider that the 
data Cfc are taken as samples from the time signal €.-'^\{t) at sam- 
pling period T where 1(0 is the unit step function, zero for negative 
<, and one for positive t. Then = e-'^'^\{kT). The z-transfonn of 
this is 



A:=— OO 0 
0 



z — e 



We wiU return to the analysis of signals and development of a table of 
useful z-transforms m Section 2.5; we first examme the use of the trans- 
form to reduce difference equations to algebraic equations and technique 
for representmg these as blodc diagrams. 

2.3.2 The IVansfer Function 

The z-transform has the same role m discrete systems that the Laplace trans- 
form has m analj^is of contmuous systems. For example, the z-transforms 
for and ti* m the difference equation (2.2) or in the trapezoid mtegration 
(2.7) are related m a simple way that permits the rapid solution of Imear, 
constant difference equations of this kind. To find the relation, we proceed 
hy direct substitution. We take the definition given by (2.8) and, in the same 
way, we define the z-transform of the sequence {ufc} as 

Uiz)= 53 ttfcz-* (2.9) 

*=-oo 
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Now we multiply (2.7) by z"* and sum over k. We get 

£ u,z''^ f: ti.-i^-* + |f E e,z-'+ f: Ciz-M. (2.10) 

kZ^oo \fe=-oo *=-oo / 

From (2.9), we recognize the left-hand side as U{z). In the first term on the 
right, we let fc 1 = j to obtain 

f: ut^iz-' = f; ujz-^^'^ = z-'Uiz). (2.11) 

k=-oo i=-oo 

By similar operations on the third and fourth terms we can reduce (2.10) to 

U{z) = z''U{z) + ^[E{z) + z'^E{z)], (2.12) 

Equation (2.12) is now simply an algebraic equation in z and the functions 
U and E. Solving it we obtain 



(2.13) 



We define the ratio of the transform of the output to the transform of the 
input as the transfer function, H{z). Thus, in this case, the transfer function 
for trapezoid-rule integration is 



(2.14) 



For the more general relation given by (2.2), it is readily verified by the same 
techniques that 



m ^ bo + biz-^ + "- + b,nZ~"' 
and if n > m, we can write this as a ratio of polynomials in z as 



+ aiz«-i + a2«''-2 J- ... -I. /r 
a{zy 



(2.15) 



22 CHAPTER 2 SYSTEMS ANALYSIS 













C^(z) = z-"£(z)* 



Figure 2.3 The unit delay. 

The general input-output relation between transforms with Unear, constant 
difference equations is ' 



U{z) = H{z)E(z). 



(2.16) 



Although we have developed the transfer function with the z-transform, it 
is also true that the transfer function is the ratio of the output to the input 
when both vary as z*. 

Because H(z) is a rational function of a complex variable, we use the 
terminology of that subject. Suppose we caU the numerator polynomial b(z) 
and the denominator a(z). The places in z where b{z) = 0 are zens of the 
transfer function, and the places in z where o(z) = 0 are the poles of H(z) 
If zo is a pole and (z - zoyH{z) has neither pole nor zero at zo, we say that 
H{z) has a pole of order p at zo. If p = 1, the pole is simple. The transfer 
function (2.14) has a simple pole at z = 1 and a simple zero at z = -1. 

We can now give a physical meaning to the variable z. Suppose we let 
all coefficients in (2.15) be zero except 6i and we take 6i to be 1. Then 
H{z) = z .But H{z) represents the transform of (2.2), and with these 
coefficient values the difference equation reduces to 



Uk = Ck-l. 



(2.17) 



The present value of the output, Uk, equals the input delayed by one period. 
Thus we see that a transfer function of z"! is a delay of one tune unit 
We can picture the situation as in Fig. 2.3, where both tune and transform 
relations are shown. 

Since the relations of (2.7), (2.14), and (2.15) are all composed of delays, 
they can be expressed in terms of z-M Consider (2.7). In Fig. 2.4 we illustrate 
the difference equation (2.7) using the transfer fonction z"* as the symbol 
for a unit delay. 

We can follow the operations of the discrete integrator by tracing the 
signals through Fig. 2.4. For example, the present value of is passed to 
the first summer, where it is added to the previous value ek-u and the sum 
IS multiplied by T/2 to compute the area of the trapezoid between Ck-i and 
e*. This IS the signal marked a* in Fig. 2.4. After this, there is another sum, 
where the previous output, Uk-i, is added to the new area to form the next 
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Figure 2.4 A block diagram of trapezoid integration as represented by (2.7). 

value of the integral estimate, u,. The discrete mtegration occurs m the loop 
with one delay, z^, and unity gain. ^ 

2.3.3 Block Diagrams and State- Variable Descriptions 
Because (2^16) is a linear algebraic relationship, a system of such relations 
f^"^^ by a system of Imear equations. These can be solved by Z 

?f^tlol dT "'"''^t " ^^P^^^' '"^^^-^ of block Sa^a^ 

Sa^^i^t .t^^ manipulate these discrete-transfer-fiiction 

relationships, there are only four primitive cases: 

onhe forward path divided by one minus the loop transfer fimction 

The transfer function of an arbitrary multipath diagram is given by 
combmations of these cases. Mason's rule* can also be^ ^ 

fiinJ^fon t irrS T ^'?^*^on of (2.2). we already have the transfer 
S^n^ ^ It IS inter^ting to connect this case with a block dia- 
gram t^mg ordy stmple delay forms for z in order to see several "canonical" 

diaJ?^ ""^^^ ^ ^^'^'^^ difference equation (2.2) to a block 

diagram mvolvmg z only as the delay operator. z-\ The &st one te wS 

"Mason (1956). See FVanklin. PoweU. and Emami-Naeini(1986) for a discussion. 
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Kg«re2.6 Block diagram of paraUel blocks. 



where 



and thus 



f ■ E{z) 
^ o(z) 



a(z)f = E(z). 



(2.18) 
(2.19) 
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KgBre2.6 Block diagram of cascade blocks. 
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Figure 2.7 Feedback transfer function. 

We can write (2.18) as 

= e - oiz^^ - a2zi - 
i{k + 3) = e{k) - aie(fc + 2) - a2^{k + 1) - asCCik). (2.20) 

Now assume we have z^^, which is to say that we have ^{k + 3) because 
is an advance operator of three steps. If we operate on this with z"^ three 
times in a row, we will get back to ^(fc), as shown in Fig. 2.8(a). Prom (2.20), 
we can now compute z^^ from e and the lower powers of z and ^ given in the 
block diagram; the picture is now as given in Fig. 2.8(b). To complete the 
representation of (2.18) and (2.19), we need only add the formation of the 
output u as a weighted sum of the variables , z^^^ ^ according to 

(2.19). The completed picture is shown m Fig. 2.8(c). 

In Fig 2.8(c), the internal variables have been named xi, X2, and X3. 
These variables comprise the state of this dynamic system in this form. 
Havmg the block diagram shown m Fig. 2.8(c), we can write down, ahnost by 
inspection, the difference equations that describe the evolution of the state, 
again using the fact that the transfer function z~^ corresponds to a one-imit 
delay. For example, we see that xz{k + 1) = X2{k) and X2{k + 1) = xi{k). 
Finally, expressmg the sum at the far left of the figure, we have 



xi(fc + 1) = -aia:i(fcj - azx'^ - + e(fc). 

We collect these three equations together in proper order, and we have 



xx{k + 1) = -oiXi(A;) - a2X2{k) - a3X3(A;) + e(fc), 
X2{k-\-l) = xx{k), 
xz{k + 1) = X2{k). 



(2.21) 
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r-1 



(c) 



SfT^ux** . diagram development of control canonical form, (a) Solving for 
i{k); (b) solvmg for f (fc + 3) from e{k) and past ^'s; (c) solving for C/(Jfc) from^'s. 
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Using vector-matrix notation/ we can write this in the compact form 
x(A: + l) = Acx(A;) + Bee(ib), 



where 



X = 




oi -02 -03 1 
10 0, 
5 1 0 J 



and 



(2.22a) 



(2.22b) 



catch a// paths by which the state variables combine in the output. The 
prob em .s caused by the bo term. If 6„ = 0. then u = b^., + to + 63x3 
and the correspondmg matrix form is immediate. However, if bo is not 0 
for example not only reaches the output through b^ but abo by the pamUe 
path with gam -boai. The complete equation is 

t* = (61 - Oi6o)xi + (62 - a2bo)x2 + (63 - 0360)2:3) + boe. 
In vector/matrix notation, we have 



u = CcX + DcC, 



where 



Cc-[6i-ai6o 62-0260 63-0360), (2.23a) 
^c = [*oJ. (2.23b) 

LVl^* ^'l' ^f"^"' ^ f of matrices. The results we require and 

^^f^ '"^'^S '"^ C. To distinguish vSors and. 

matrices from scalar variables, we will use bold-face type. 
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We can combine the equations for the state evolution and the output to give 
the very useful and most compact equations for the dynamic system, 



x(A; + l) = Acx(A:) + Bce(fc), 



(2.24) 



where Ac and Be for this control canonical form are given by (2.22), and 
Cc and Dc are given by (2.23). 

The other canonical form we want to illustrate is called the observer 
canonical form and is found by startmg with the difference equations in 
operator/ transform form as 

z\ + aiz^u + a2zu + a^u = bozh + biz^e + b2ze + 63c. 

In this equation, the external input is e(A:), and the response is li(ifc), which 
is the solution of this equation. The terms with factors of z are tune-shifted 
toward the future with respect to k and must be elimmated in some way. 
To do this, we assume at the start that we have the u{k), and of course the 
e(A:), and we rewrite the equation as 

he - a^u = z\ + aiz\ + a2zu ~ h^z^e - biz^e - 62^6. 

Here, every term on the right is multiplied by at least one power of z, and 
thus we can operate on the lot by z'^ as shown in the partial block diagram 
drawn m Fig. 2.9(a). 

Now m this internal result there appear a2U and -626, which can be 
cancelled by adding proper multiples of u and e, as shown in Fig. 2.9(b), 
and once they have been removed, the remainder can again be operated on 
by z"^. 

li we continue this process of subtracting out the terms at * and oper- 
ating on the rest by z'^, we finally arrive at the place where all that is left 
is u alone! But that is just what we assumed we had in the first place, so 
connecting this term bade to the start finishes the block diagram, which is 
drawn in Fig. 2.9(c). 

A preferred choice of numbering for the state components is also shown 
in the figure. FoUowing the technique used for the control form, we find that 
the matrix equations are given by 



x(A; + 1) = Aox(fc) + Boe(fc), 
«(A:) = Cox(*:)+Doe(A:), 



(2.25) 
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zu — b^e — by€ 



(b) 




(2.25) 



Figure 2.9 Block diagram development of observer canonical form, (a) The first 
partial sum and delay; (b) the second partial sum and delay; (c) the completion 
with the solution for u{k). 
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so 
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z^- 0.5H"Q.2S 
z*- 1.62-^0.81 



Figure 240 Block diagram of a cascade realization. 



where 



Ao = 

Bo = 

Co 
Do 



-ai 1 0 
-02 0 1 
-as 0 0 

6i — 6ooi 

63 — 6003 . 
[1 0 0], 
[hoi 



The block diagrams of Figs. 2.8 and 2.9 are called direct canonical re- 
alizations of the transfer function H{z) because the gains of the realizations 
are coefficients in the transfer*function polynoniials. Another useful form is 
obtained if we realize a transfer function by placing several first- or second* 
order direct forms in series with each other, a cascade canonical form. In 
this case, the H{z) is represented as a product of factors, and the poles and 
zeros of the transfer function are clearly represented in the coefficients. 

For example, suppose we have a transfer function 

. + 0,5z^ - 0-25Z + 0.25 
z4.2.6z3 + 2.4z2-0.8z 
^ + 0.5^ + 0.25) 

(z2^z)(z2-1.6z + 0.8)' 

The zero factor z+1 can be associated with the pole factor z^— z to form one 

seeoM^^rda^^tem^ aimt'^^zm^^f^ 

with the second-order pole factor z^-1.6z+0.8 to form another. The cascade 
factors, which could be realized in a direct form such as control or observer 
form, make a cascade form as shown in Fig. 2.10. 

2.3.4 Relation of Transfer Function to Pulse Response 

We have shown that a transfer function of z""^ is a unit delay in the time do- 
main. We can also give a time-domain meaning to an arbitrary transfer func- 
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Table 2.1 Step-by-step construction of the unit 
pulse response for Fig. 2.4. 



k 












0 


0 


1 


r/2 


0 


T/2 


1 


1 


0 




T/2 


T 


2 


0 


0 


0 


T 


T 


3 


0 


0 


0 


T 


T 



tion. Recall that the z-transform is defined by (2.8) to be E{z) = Eefcz"* 
and the transfer function is defined from (2.16) as H{z) when the input and 
output are related by U{z) = H{z)E{z). Now suppose we deUberately select 
e{k) to be the unit discrete pulse defined by 



n. (* = 0), 
^* lo, (fc#0), 



Then it follows that E{z) = 1 and therefore that 

U{z) = H{z). 



(2.26) 



(2.27) 



Thus the transfer function H{z) is seen to be the transform of the response to 
a umt-pulse mput. For example, let us look at the system of Fig. 2.4 and put 
a unit pulse in at the Cit-node (with no signals in the system beforehand)." 
We can readUy foUow the pulse through the block and build Table 2.1. 

Thus the unit-pulse response is zero for negative k, is T/2 at = 0, and 
equals T thereafter. The z-transform of this sequence is 

H(z) = f;«*z-*^f;fe,^-* 



In this development we assume that (2.7) is mtended to be used as a formula for 
computing values of «ik as t increases. There is no reason why we could not also 
solve for u* as * takes on negative values. The direction of time comes from the 
application and not from the recurrence equation. 
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+63/10^"^ 

+62^1 
+61^2^"^ 



coefficients, 



Flgnre2.ll Representation ofthe product £?(z)^(z) as aproduct of polynomials. The extrapol 



If we add T/2 to the ^o.tenn and subtract T/2 from the whole series 
we have a simpler sum, as follows: ' 



ik=0 



2 



_ 2r-r(i-z-i) 
2(1 -z-i) 

2(1 -z-l) 
2z-l 



(1 < \z\) 



(1 < kl). 



(2.28) 



By extension, 
be +00: 



Negative valui 
equals zero. \ 
is nonzero for 
spends before t 
discrete convol 
relates input a 
^tems. 

Tb verify (: 



anai^is^r^hi^^i:^:^ ^'-'^^ ^-^^ 

fanction IS obt^ed by multii^lying the mfinite polynomiSnf ^^^^ 
fh^tZ uT'f " ''''' m4r^ion,te wiU ^i^e 

1 7hP ioi f^"?^ "^^^^ P^^i that i^rodue^^-S 

IS the polynomial product given m Fig. 2.11. P^auces u (z) 

Smce this product has been shown to be V(z) = if m„cf *k 

fore foUow that the coefficient of in JplductTu^lS^ 



Nowlet fc-j = 
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coefficieiits, we have the relations 

^l = eohi + eiho, 

U2 = €0^2 + eiftl + 62/10, 

The extrapolation of this simple pattern gives the result 

k 

j=o 

By extension, we let the lower limit of the sum be —00 and the upper limit 
be +00: 

00 

ti/k= E ^i'^^-i- (2.29) 

i=-oo 

Negative values of j in the siun correspond to inputs applied before time 
equals zero. Values for j greater than k occur if the tmit-pulse response 
is nonzero for negative arguments. By definition, such a system, which re- 
sponds before the input that causes it occurs, is called noncausaL This is the 
discrete convolution sum and is the analog of the convolution integral that 
relates input and ipipulse response to output in linear, constant, continuous 
systems. 

To verify (2.29) we can take the z-transform of both sides: 

00 00 00 

jk=-oo jk=-oo j=-oo 

Interchan^ng the smn on j with the sum on k leads to 

3— —00 k=—oo 

Now let ib — J = I in the second sum: 
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but = z'^z'^y which leads to 

j=—oo /=— OO 

and we recognize these two separate sums as 

U{z) = E{z)Hizy 

We can also derive the convolution sum from the properties of linear- 
ity and stationarity. First we need more formal definitions of "linear'' and 
"stationary." 

1. Linearity: A system with input e and output tx is linear if superposition 
applies, which is to say, if ui{k) is the response to ei{k) and U2{k) is 
the response to 62 (fc), then the system is linear if and only if, for every 
scalar a and ^, the response to aei + /?e2 is aui + fiu2. 

2. Stationarity: A system is stationary ^ or time invariant, if a time shifk in 
the input results in only a time shift in the output. For example, if we 
take the system at rest (no internal energy in the s)rstem) and apply a 
certain signal e(A;), suppose we observe a response u{k). K we repeat this 
experiment at any later time when the system is again at rest and we 
apply the shifted input, e{k — TV), if we see u{k — iV), then the system 
is stationary. 

These properties can be used to derive the convolution in (2.29) as fol- 
lows. If response to a imit pulse at A: = 0 is h{k)j then response to a pulse 
of intensity eo is eoh{k) if the system is linear. Rirthermore, if the system is 
constant, then a delay of the input will delay the response. Thus, if 



(eh 
0, 



then the response will be eihk^i. 

Finally, by ffieMfy again, the totsd response at time i to a sequence of 
these pulses is the sum of the responses, namely, 

Ufc = eohk + eihk^i H h eihk-i + h Cfe/io> 

or 
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Now note that if the input sequence began in the distant past, we must 
mclude temas for / < 0, perhaps back to / = -oo. Similarly, if the system 
should be noncausal, future values of c where / > k may also come in. The 
general case is thus (again) 

oo 

«*= 53 ^if^k-t- (2.30) 

f=-oo 

2.3.5 External Stability and Jury's Itest 

A very important qualitative property of a dynamic system is stability, and 
we can consider internal or external stabiUty. Internal stabUity is concerned 
with the responses at aU the internal variables such as those that appear at 
the delay elements in a canonical block diagram as in Fig. 2.8 or Fig. 2.9 (the 
state). Otherwise we can be satisfied to consider only the external stability 
as given by the study of the input^utput relation described for the Unear 
stationary case by the convolution (2.30). These differ m that some internal 
modes might not be connected to both the input and the output of a given 
system. 

For external stability, the most common definition of appropriate re- 
sponse is that for every Bounded Input, we should have a Bounded Output. 
If this is true we say the system is BIBO stable. A test for BIBO stability can 
be given directly in terms of the unit-pulse response, fc*. First we consider 
a sufficient condition. Suppose the input ek is bounded, that is, there is an 
M such that 

|ej|<M<oo for all/. (2.31) 

If we consider the magnitude of the response given by (2.30), it is easy to 
see that 

\uk\<\Y.eihk-i\, 
which is surely less than the sum of the magnitude as given by 

<fl\ei\\hk-i\. 

— OO 

But, because we assume (2.31), this result is in turn bounded by 

oo 

<MY,\hk-l\. (2.32) 
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Thus the output will be bounded for every bounded input if 

oo 

£ \hk-i\ <oo. (2.33) 

(=-oo 

This condition is also neceissary, for if we consider the bounded (by 1!) input 

= 0 {h-i = 0) 

and apply it to (2.30), the output at fc = 0 is 



tio= 53 

/=-00 



= f; \h-i\. (2.34) 

l=-oo 

Thus, unless the conditibn given by (2.34) is true, the system is not BIBO 
stable. 

Example 2.3: The test given by (2.34) can be applied to the unit 
puke response used to compute (2.13) and ^ven as the Ufc-colunm 
in Table 2.1 on page 31: 

ho = T/2, 

hk = T, k>0, 

oo 

£ \hk\ = r/2 + J3 T = unbounded. (2.35) 
1 

Thus this discrete approximation to integration is not (BIBO) stable! 
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f 

(2.33) 
ed (by 1!) input 



(2.34) 
;m is not BffiO 



plied to the unit 
) the Ufc-column 



(2.35) 
t (BmO) stable! 



Example 2.4: As a second example, we consider the difference 
equation (2.2) with aU coefficients except oi and 6o equal to zero: 

ttfc = oittfc_i + boek. (2-36) 

The unit-pulse response is easily developed from the first few terms 
to be 



uk = hk = boa!', > 0. 
Apptying the test, we have 



(2.37) 



— oo 00=0 



=: unbounded 



(|a| < 1) 
(|a| > !)• 



Thus we conclude that the system described by this equation is BffiO 
stable if |o| < 1, and unstable otherwise. 

For a more general rational transfer function with many simple pol^, 
we c^ expand the function in partial fractions about its poles and the 
Zr^nLg pulse response wiU be a smn of respective terms As we saw 
eL^^a^le is insile the unit circle, the corresponding P^,^^^^^ 
dLys' withlime geometricaUy and is stable T^us, if all po^ axe ms^d th 
unit circle the system with rational transfer function is stable, if at least 
Te po" Von of outside the unit circle, the -rrespondi.^ sj^tem ^^^^^ 
BffiO stable. With modern computer programs avai able, findmg the pol^ 
nf . Lticular transfer function is no b ig deal. Sometimes, however, we w^h 

system, Jbe^gji 

SLuois case, such'a test was provided by R^uth; in the <ii--te^^^^^^^^^ 
most convenient such test was worked out by Jury Blanc^rd(1961^^ 
The Jury test is in the same spirit as the Routh test (see FVanWm, 
Powell andEmami-Naeini for a discussion) in that we form two rows from 
lir^o^ents of length n. and from these, by a seri^ --^^'^^^^ . 
terminants, we compute a successor row of length n - 1. With this reduced 



I 
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length row we form another successor of length n - 2 and so on until we 
have a row of length 1, The test consists of examining the sign of the first 
entries in selected rows. As with the Routh test, the Jury test is much more 
difficult to derive than to use; here we illustrate only the use. 

If we have a transfer function H{z) = b(z)/a{z), then this system will 
be stable if and only if all roots of a(z) = ooz" + aiz"""^ + • • • + a„ are inside 
the unit circle- To test for this condition by the Jury test, multiply a{z) 
by -1 if necessary to make the sign of ao positive. Then form rows of the 
coefficients, the even rows being in reversed order, as follows: 



Ex 

the 



Tb 



Oo ax 



On 



The entries in the third row are formed from the second-order deter- 
mmants using the first column of the first two rows with each of the other 
columns from these rows starting from the right and dividing by oq. The 
result can be expressed by the formulas: 



froi 



Rro 



. "n 
60 = OO On, 

I. ^ 

01 = Oi On-l, 

The elements in the thfrd row are reversed to form the fourth row and the 
process is repeated. For example, the elements of the fifth row are given by 

The original polynomial is stable ( has all roots inside the unit circle) 
if all the terms in the first columns of the odd rows are positive, that is, if 
ao>0,6o>0,co>0,.... 

This test is readily implemented in a computer program.^ 



and 



FVo] 
Fig 



^See STABLE m Table E.l in Appendix E. 



Figure 2.1 
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a{z)^z^ + aiz + a2. 



The Jury array is 



<»1 02 

oi 1 



1 

l-a| 
ai — ai02 

(l-ai)^-o?(l-a,)2 

Rom r<w three we have the condition that 1 - > o 
from this we conclude that ' 

-1 < 02 < 1. 

FVom row five, we can factor out (1 - aj)^ to conclude that 
(l + 02)2>a,2 

and thus, 

02 + 1 > ai and 02 + 1 > -aj. 
PVom these inequaUties, we can draw the 5*06^% triangle shown 



and 



in 




Figure 2.12 StabiUty triangle for the general, real, second-order 



polynomial. 
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Example 2.6: For the next example, consider the polynomial 
a{z) = z^- %\z^ + l.6z - 0.4. 
The Jury array is 



1 

-0.4 
0.84 
0,76 
0.1524 
-0.139 
0.0256 



-2.1 
1.6 

-1.46 
-1.46 
-0.139 
0.1524 



1.6 
-2.1 
0.76 
0.84 



-0.4 
1 



nrwlsi r ""^^ ^ > > 0,0.1524 > 0, 

m > 0, and we conclude that a system with this polynomial as 
ite denwnmator would be stable. As a matter of feet, the poles 
atz = 0.5and2 = 0.8±0.4i. 



are 



Example 2.7: As a third example, consider the polynomial 

z*-2.6ar* + 2.4?-0.8 
The odd rows only of the Jury array are 

1 -2.6 2.4 -0.8 

0.36 -0.68 0.32 . 

0.0756 -0.0756 
0 . . . 

Ftom thae computations we conclude that the polynomial does 
not napsaiytSilMsatsi^ ^^^^^^^^mm^ 
is ^ and a«naU perturbation would send it either way, mside or 
outode there mm be at least one root exactly on the Lt circle, 
b fact, the roots are z = 1 and 2 = 0.8 ± 0.4i. 



As an aid to tatmg stability, it can be shown that it is necessary for a 
stable polynonual (with positive first term) that the polyno^S^^eS 
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at 2: = 1 and z = -1 must both be positive. The first value is just the sum 
of the coefficients, and the second is the sum with alternating sign changes. 
These two tests are quickly done and can save time if the only purpose is 
to be sure that the system is stable; many imstable systems will be rejected 
by these simple tests without gomg through the entire Jury array. In the 
previous case, for instance, the sum of coefficients is zero and we need go no 
further; the polynomial cannot be stable. 



2.4 DISCRETE MODELS OF SAMPLED-DATA 
SYSTEMS 

The systems and signals we have studied thus far have been defined in 
discrete time only. Most of the dynamic systems to be controlled, however, 
are continuous systems and, if linear, are described by continuous transfer 
functions in the Laplace variable 5. The interface between the continuous 
and discrete domains are the A/D and the D/A converters as shown in 
Fig. 1.1. In this section we develop the analysis needed to compute the 
discrete transfer function between the samples that come firom the digital 
computer to the D/A converter and the samples that are picked up by the 
A/D converter. ^° The situation is drawn m Fig. 2.13. 

2.4.1 Using the zrTVansform 

We wish to find the discrete transfer function firom the input samples u{kT) 
(which probably come firom a computer of some kind) to the output samples, 
y{kT) picked up by the A/D converter. Althotigh it is possibly confusing 
at first, we follow convention and call the discrete transfer function G{z) 
when the continuous transfer function is G{s). Although G{z) and G{s) are 
entirely different functions, they do describe the same plant, and the use 
of 5 for the continuous transform and z for the discrete transform is always 
maintained. To find G{x) we need only observe that the y(fcr) are samplesnof 
the plant output when the input is firom the D/A converter. As for the D/A 
converter, we assmne that this device, commonly called a zero-order hold 
or ZOH, accepts a sample u{)iT) at t = fcr and holds its output constant 



In Chapter 3, a comprehensive fi-equency analysis of sampled data systems is 
presented. Here we undertake only the special problem of finding the sample-to- 
sample discrete transfer fimction of a continuous system between a D/A and an 
A/D. ^ 
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2.4 DISC 
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y(A:r) 



Figure 2.13 The prototype sampled-data system. 



at this value until the next sample is sent at t = kT + T. The piecewise 
constant output of the D/A is the signal, u(*), that is applied to the plant. 

Our problem is now really quite simple because we have just seen that 
the discrete transfer function is the z-transform of the samples of the output 
when the input samples are the unit pulse at fc = 0. If u{kT) = 1 for fc = 0 
and u{kT) = 0 for / 0, the output of the D/A converter is a pulse of width 
T seconds and height 1, as sketched in Fig. 2.14. Mathematically, this pulse 
is given by 1(<) - l{t - T). Let us call the particular output in response to 
the pulse shown in Fig. 2.14 yi{t). This response is the diflFerence between 
the step response [to l{t) J and the delayed step response [to l{t - T)]. The 
Laplace transform of the step response is G{s)/s. Thus in the transform 
domain the unit pulse response of the plant is 



yi(5) = (l-e-^*)^, 



(2.38) 



and the required transfer function is the ^-transform of the samples of the 
inverse of ¥1(3) , which can be expressed as 

G(z)^Z{yi(kT)} 

= ^{^-Hil(3)}} = 2{yi(^)} 

This is the sum of two parts. The fiirst is 2^{G{s)/s}, and Uiie second is 

Z{e'^^G{s)/s} = z-'Z{G{s)/s} 
because e"^« is exactly a delay of one period. Thus the transfer function is 
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(2.39) 
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Figure 2.14 D/A output for unit-puke input. 

Example 2.8: As a first example of computing such discrete trans- 
fer functions, suppose G{3) = a/{s + a). Then 

Gjs) ^ Q ^ 1 1 
s s{s + a) 3 5 + a' 

and the corresponding time function is 



£-i|^| = l(i)-e-'l(0. 



The samples of this signal are l(*:r) - e-'^l{kT), and the z- 
transfonn of these samples is 



z(l-e-°^) 



We''6oM#have»g3ne»to the4-abl&^B> J^ppeacUx^Bvsm 

directly as entry 12. Now we can compute the desired transform by 

applying (2.39) 



G(z) = 



z-1 z(l-e-«^) 
z (z-l)(z-e-«^) 




mm 
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Example 2.9: We consider the double integrator characteristic of 
a single mass, such as the sateUite, for which the transfer function is 
G{s) = We have 



This tinie we refer to the tables m Appendbc B and find that the 
transform of is 

2 (z-l)3' 

and therefore 

= 2(7ri)^- (2.41) 



For more complex systems than those in Examples 2.8 and 2.9, use of a 
CAD package is recommended.^^ 

2.4.2 Continuous Time Del^y 

We now consider computmg the discrete transfer function of a continuous 
system with pure time delay. The responses of many chemical process- 
control plants exhibit pure time delay because there is a finite time of trans- 
port of flmds or materials between the process and the controls and/or the 
sensors. Also, we must often consider finite computation time m the digital 
controUer, and this is exactly the same as if the process had a pure time 
delay. With the techniques we have developed here, it is possible to obtain 
musteat'^^*^ *^^*er fimction ^^^uch prwesses exactly, as Example 2.10 

Example 2.10: We consider the example suggested by the fluid 
nuxer problem described m Appendbc A.3, for which 



G(s) = e-^His) 



"See X-C2D in Table E.l. 
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The term e~^^ represents the delay of A seconds, which mcludes both 
the process delay and the computation delay, if any. We assume that 
H{3) is a rational transfer function. To prepare this function for 
computation of the 2-transform, we first define an integer £ and a 
positive number m less than 1.0 such that A = f T - mT. With these 
definitions we can write 

— — — c . 

5 S 

Because ^ is an integer, this term reduces to z"^ when we take the z- 
transform. Because m < 1, the transform of the other term is quite 
dnect. We select H{s) = a/{s + a) and, after the partial firaction 
expansion of H{s)i we have 

To complete the transfer function, we need the z-transforms of the 
inverses of the terms in the braces. The first term is a unit step 
shifted left by mT seconds, and the second term is an exponential 
shifted left by the same amoimt. Because m < 1, these shifts are less 
than one ftdl period, and no sample is picked up in negative time. 
The signals are sketched in Fig. 2.15. 

The samples are given by l{kT) and e-^W'»)i(fcr). The cor- 
respondmg z-transforms are z/(z - 1) and ze'^^'^^/Cz - e"-**^). Con- 
sequently the final transfer function is 

^ z-l ( z[z- e-"^ -{z- 1)6-°"*^] \ 

^z ^ l^^z ^e^'^^^ f ■ 

where the zero position is at -a = -(e-""*^ - e-«^)/(l - cy*"*^). 
Notice that this zero is near the origin of the z-plane when m is near 
1 and moves outside the unit circle to near -oo when m approaches 
0. For the specific values of the mixer, we take o = 1, T = 1, and 
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Figure 2*15 Sketch of the shifted signals showing sample points. 

A = 1.5. Then we can compute that < = 2 and m = 0.5. For these 
values, we get 



G{z) = 



z + 0.6065 
z2(z^ 0.3679)' 



(2.42) 
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2.4.3 State-Space Form 

Computing the z-transform using the Laplace transform as in (2.39) is a very 
tedious business that is unnecessary with the availability of computers. We 
will next develop a formula tising state descriptions that will remove most 
of the calculations to the computer, where it is better done. A continuous, 
linear, coxistant-coefficient system of dififerential equations can always be 



x = Fx + Gti + Giw, 



(2.43) 



where u is the control input to the system and w is a disturbance input. 
The output can be expressed as a linear combination of the state, x, and 
the input as 



The represc 
representation, 
( = Tx is also J 

If welet( = 



If we designate 
then 



y = Hx + Ju. 



(2.44) 
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Figure 2.16 Satellite attitude control in classical representation. 

Often the sampled-data system being described is the plant of a control 
problem, and the parameter J in (2.44) is zero and will frequently be omit- 
ted. 



Example 2*11: Application of state representation to the equa- 
tions of the satellite attitude-control example shown in Fig. 2.16 and 
described in Appendix A yields 



P G 



(2.46) 



which, in this case, turns out to be a rather involved way of writing 

0 = u. 

The representations (2.43) and (2.44) are not unique. Given one state 
representation, ai^ nonsingular linear transformation of that state such as 
{ = Tx is also an allowable alternative realization of the same system. 

If we let ^ = Tx in (2.43) and (2.44), we find 

C = Tx = T(rx + Gu + Gittf) 

^ = TFT-*C + TGu + TGiu;, 
y = HT-^e+«'«- 

If we designate the system matrices for the new state ^ as A, B, C,andD, 
then 

^= A^ + Bu + Biw, y = Ce + Dti, 
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^^2.17 System definition with «.«™ I- 

"on with samphng operations shown. 



where 



i^^,?;??^ 2.12: As an illustration,^ can let ^ 
in (2.45); or, in matrix notation f k1 T ^ = ^2 and = x, 
the states is "oiauon, the transformaf{«», ^ . 



transformation to interchange 



In this case T"^ = t v 



gives 

C = fO IJ. 



equat on could be represented by^^ItS^ dl^''^^ '^^^ ^'S^-°^<»«- differ?^ 
ca^omcal fonn. Aiso, there is ZT^m^T^ " ^"^^observ^ 
the partial-fraction expansion ofTt?!^ a ^''"P"°"^°"e«Pondingto 
-n take a general t^^^Trt ct^^^^"" '-nsLmatioiL 

and, subject to some technical restricHont ''^^^^^^^^^ « discrete system 
WetLT ^orms,rSr' ^ Won in 

obtainmgthe ^fciS^^atio^^^^ - ^-eral method for,^ 

uoos plant. Fig. 2.17 ag^ deSJ^e 3 *5" t^e contin- 

eration Ulthnately, the^X^ntr^er^Taf .7 ^^*""^ 
on that sequence by means of a dSe^^ . ^ ' '^"^^^^ ^^W. operate 
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We wiU solve the general equation in two steps. We begin by solving 
the equation with only initial conditions and no external input. This is the 
homogeneous equation 



Xh = Fxa(<), Xhito) = xo 



(2.46) 



To solve this, we assume the solution is sufficiently smooth thiatt a series 
expansion of the solution is possible: 



Xh{t) = Ao + Ai{t - <o) + A2{t - hf + ' 



(2.47) 



if we let t = fo, we find inmiediately that i4o = xq. If we differentiate (2.47) 
and substitute into (2.46), we have 

ill + 2i42(< - <o) + 3A3(< - <o)^ + • • = Fx/i 

and, at t = to, Ax = Fxq. Now we continue to differentiate the series and 
the differential equation and equate them at to to arrive at the series 



I + F(t-to) 



2 6 
This series is defined as the matrix exponential and written 

Xfc(<) = e^(*-*o)x(io), 
where, by definition, the matrix exponential is 

eF(t-to) = I + F(t - to) + f2 + + 



Xq. 



(2.48) 



oo 



k=0 



k\ • 



(2.49) 



It C2QiWitt6wiItMt''ttt^ §S]t3:t%tt^veB (2.48) is uiiiqute, wMmfe^«to 
very interesting properties of the matrix exponential. For example, consider 
two values of f : t\ and t2- We have 



x(ti) = c^t*i-*o>x(to) 



and 



X(t2) = C^<*2-*»>x(«o). 




ion 
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Because <o is arbitrary also, we can express xffo) as iff!,. _ 
began at ti, for which ^ " equation soluti, 

Substituting for x(ti) gives 



(2.50) 



C*'(*2-<0) _ eF«2-<l)gP(«l-<o) 

for aU <2, <i, to. Note especiaUy that if t2 = to, then 

I = c-*'(*i-«o)cP(«i-«o) 
Thus we can obtain the mverse of e*"* hv .r,., i u . 

we obtain •~"™«<™ *W m (2.48)J. Sub««tultag (2.51) into (2.43), 

Pe'C-Wy + ^ Fe»('-'.)v + G„, 

and, usmg the feet that the inverse is fo.,nri k. ». • 
exponent, we can solve for v ^ "^^^^ the 

v(i) = e-*'(*-^)Gtt(e). 



r="rbJ*SS*^!f^«'*°««'^«'«* mathematician a?^ w 
w - u, but because the equations are lln*.ar «. t . , t^'^o-l8l3). We assume 
«I auons ate hnear, the effect of «; can be added later. 



The total solu 
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Assuming that the control u{t) is zero for * < to> we can integrate v from to 
to t to obtain 

Hence, from (2.51), we get 

and simplifying, using the results of (2.50), we obtmn the particular solution 
(convolution) 

Xp(<) = /' e^(*-'^)Gti(r)rfr. (2,52) 
•'to 

The total solution for «; = 0 and ix # 0 is the sum of (2.48) and (2.52): 

x(t) = e^(*-^o)x(<o) + f e^(*-^)Gu(r)dr. (2.53) 

We wish to use this solution over one sample period to obtain a difference 
equation: hence we juggle the notation a bit (let t-kT + T and to equal 
kT) and arrive at a particular version of (2.53): 

AkT + T) = e^x{kT) + e^(«=^+^--)Gu(r)(ir. (2.54) 

This result is not dependent on the type of hold because u is specified 
in terms of its continuous time history, u{t), over the sample mterval. A 
common and typically vaUd jij^jamptiQn is that of a zecohOEder hold (ZOH) 
with no delay, that is, 

tt(r) = tt(fcT), kT<T<kT + T. 

If some other hold is implemented or if there is a delay between the appli- 
cation of the control from the ZOH and the sample point, this fact can be 
accounted for in the evaluation of the integral m (2.54). The equations for a 
delayed ZOH will be given in the next subsection. To facilitate the solution 
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of (2.54) for a ZOH with no delay, we change variables in the integral from 
Ttofj such that ^ 



fl=zkT + T-r. 



Then we have 



If we define 



zikT + r) = e^^^ikT) + ]\^HnGu{kT). (2.55) 



(2.56a) 

T = ]\^nir,G, (2.56b) 
Eqs. (2.55) and (2.44) reduce to difference equations In standard form: 
x(A: + 1) = #x(A;) + rtt(ib) + ri«;(ib), 



y(*) = Hx(ifc), 



(2.57) 



where we include the effect of an impulsive or piecewise constant disturbance 
K ' f^ltff ./ = 0 in this case. If «; is a constant, then T, is rive.; 
by (2.56b) with G replaced by Gi. ti; is an impulse, then T, = G, " The 
9 series expansion x i- xro 



can also be written 



whiere 



2! 



# = I + FT*, 



+ • • • 

3! ^ ' 



(2.58) 



If «;(<) is not a function of only its sample values, then an mtegral like that of 



an 



1. 
2. 
3. 
4. 
5. 

6. 
7. 
8. 
9. 
10. 

Figure 2. 
(The left 

The r ii 



We( 



which hi 
then fint . 
of 7? an< 
Kallsti^j t 
by Mole 
andrfo 
we know 
F, G,ar 

"See X-( 
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! integral from 



1. 
2. 
3. 
4. 
5. 

6. 
7. 
8. 
9. 



Select sampling period T and description matrices F and G. 
Matrix I <- Identity 
Matrix 9 *-I 

jfc i- 11 [We are using iV = 11 in (2.60).] 

If A; = 1, go to step 9. 

FT 

Matrix * 1 + — * 
k 

k*-k-l 
Go to step 5. 
Matrix T *- r*g 



10. Matrix*<-I + Fr* 



Figure 2.18 Program logic to compute $ and T from F, G, and T for simple cases. 
(The left arrow, ♦-, is to be read "is replaced by.") 

The r integrd in (2.56) can be evaluated term by term to give 



We evaluate 9 a series in the form 



(2.59) 



, , FT/ FT/ Fr /, FT\\ \ 

which has better numerical properties than the direct series of powers. We 
then find T from (2.59) and * from (2.58). A discussion of the selection 
of N and a technique to compute * for comparatively large T is given by 
Kallstrom (1973), and a review of various methods is found in a classic paper 
by Moler and Van Loan (1978). The program logic for computation of ^ 
and r for simple cases is given in Fig. 2.18. All control design packages that 
we know of contain logic to compute ^ and T from the continuous matrices 
F, G, and the sample period T}^ 



^^See X-C2D in Table E.l m Appendix E. 
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lb compare this method of representing the plant with the discrete trans- 
fer functions, we can talte the z-transform of (2.57) with w = 0 and obtain 



therefore 



[zl-9']X{z)=^TU{z), 
n^) = HX(z); 



^=H[zI-*]->r. 



(2.61a) 
(2.61b) 



(2.62) 



Ebrample 2.13: Fbr the sateUite attitude-control example, the * 
and r matrices are easy to calculate using (2.58) and (2.59) and the 
values for F and G defined in (2.45). Since = o in this case we 
have ' 

hence, using (2.61), we obtam 
2 

which is the same result that would be obtained using (2.39) and the 
z-transform tables, »v ^ 



n.inw VI?'t*^ IT^^.*^ denominator is the deter- 

mmant det(2rl - $), which comes from the matrix inverse in (2.62). This 



determ 
the zei 
at z = 
equatic 
We 
space c 
interpn 
spondii 
tion ha 
this im] 



have a r 
for this 



which is 

nimierica 
tine.^5 

Aloni 
the systei 
Thus if ^ 
that Y{z^ 
of (2.57), 



Once moi 
the deteri 



^^See EIGI 
'^We do nc 



r 
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determinant is the characteristic polynomial of the transfer function, and 
the zeros of the determinant are the poles of the plant. We have two poles 
at z = 1 in this case, corresponding to the two integrations in this plant's 
equations of motion. 

We can explore further the question of poles and zeros and the state- 
space description by considering again the transform equations (2.61). An 
interpretation of transfer-function poles from the perspective of the corre- 
sponding difference equation is that a pole is a value of z such that the equar 
tion has a nontrivial solution when the forcing input is zero. FVom (2.61a), 
this implies that the linear eigenvalue equations 

[zI-*]X(z) = [0] 

have a nontrivial solution. From matrix algebra the well-known requirement 
for this is that det(zl — $) = 0. In the present case, we have 

<i.,.I-„ = <.e.[[j 

= (z-l)2 = 0, 



which is the characteristic equation, as we have iseen. To compute the poles 
numerically when the matrices are given, one would use an eigenvalue rou- 
tine.^5 

Along the same line of reasoning, a system zero is a value of z such that 
the system output is zero even with a nonzero state-and-input combination. 
Thus if we are able to find a nontrivial solution for X(2;o) and t/(zo) such 
that y(2io) is zero, then zo is * zero of the system. Combining the two parts 
of (2.57), we must satisfy the requirement 



(2.63) 



Once more the condition for the existence of nontrivial solutions is that 
the determinant of the square coefficient system matrix be zero.^^ For the 



tor is the deter- 
in (2.62). This 



^*See EIGENV in T^ble E.l. 

^®We do not consider here the case of different numbers of inputs and outputs. 
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satellite example, we have 



det 



2-1 

0 
1 



-T 
z-l 
0 



-r2/2 

-T 
0 



[z-l -T \ 



1) 



2^2 y2 

= + ^(z + l). 

zeros, good algorithms exist in matrix ^ebrl '^ ' *^^«^ion 
2.4.4 State-Space Models for Systems with Delay 

x(0 = Px(0 + Gtt(<-A), 
» = Hx 

The general solution to (2.64) is given by (2.54); it is 

x«) = e^(*-'o)x(«o) + f\n*-r)Gu(r - X)dr. 
Jto 

If we let /o = fcr and t = ifcr + T, then 

x(*r + T) = e"'x(*7') + c-(*^+^-^)G„(r - A) dr. 



(2.64) 



^J^Hy Khe\!^^^L^^^ ^^''f ^'^t"' ^ *o -count 



r 
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the transfer 
ransmission 



aodels from 
r present the 
2 prediction 
1 as defined 
lym control 



(2.64) 



ir. 



i to account 
hem as very 



u{t) 



-L 



kT-HT kT-iT-^T kT-n + 2T t 
Figure 2.19 Sketch of a piecewise input and time axis for a system with time delay. 

If we substitute = fcT + < — r for r in the integral, we find a modification 
of (2.55): 

x(A;r + T) = e'^x(ibr) + T e^'^GuikT + T-X- v){'dfi) 

Jt 

= e^x(ifer) + r e^^Gtt(ikr + T - A - T/) dT/. 
Jo 

If we now separate the system delay A into an integral number pf sampling 
periods plus a fraction, we can define an integer £ and a positive nmnber m 
less than one such that 



and 



X-tr-mT, (2.65) 
0 < m < 1. 

With this substitution, we find that the discrete system is described by 

x(*r + r) = e^x(Jbr)+ / e^'^GuikT + T-£T + mT'-ri)dri, (2.66) 

If we sketch a segment of the time axis near t = kT - £T (Fig. 2.19), the 
nature of the integral in (2.66) with respect to the variable tj will become 
clear. The integral runs for Tj from 0 to T, which corresponds to t from 
kT-er+T + mT backward to kT-lT + mT. Over this period, the control, 
which we assimie is piecewise constant, takes on first the value u{kT—tT+T) 



I 
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and then the value u{kT-tr). Therefore, we can break the mtegral in (2.66) 



into two parts as follows: 





= *x(fcr) + TmikT - iT) + T2u(kT --eT + T). (2.67) 



In (2.67) we defined 




To complete our analysis it is necessary to express (2.67) in standard state- 
space form. To do this we must consider separately the cases of ^ = 0, € = 1, 
and e>l. 

For ^ = 0, A = -mr according to (2.65), which implies not delay but 
prediction. Because mT is restricted to be less than T, however, the output 
will not show a sample before A = 0, and the discrete system will be causal. 
The result is that the discrete system computed with ^ = 0, m 7^ 0 will 
show the response at t = 0, which the same system with € = 0, m = 0 would 
show at < = mT. In other words, by takmg ^ = 0 and m # 0 we pick up 
the response values between the normal sampUng instants. In z-transform 
theory, the transform of the system with ^ = 0, m 7^ 0 is called the modified 
Z'transform}^ The state-variable form requires that we evaluate the integrals 
in (2.68). To do so we first convert Ti to a form similar to the integral for 

From (2.68) we factor out the constant matrbc G to obtain 




T 



H we set a = 1/ — mf in thas^ inte^ail> we have 




(2.69) 



"See Jury (1964) or Ogata (1987). 
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in (2.66) 



For notational purposes we will define, for any positive nonzero scalar num- 
ber, a, the two matrices 



(2.67) 



(2.68) 

d state- 
),^=1, 

Jay but 
output 
causal. 
= 0 wiU 
) would 
pick up 
msform 
nodified 
itegrals 
?-al for 



(2.69) 



$(a) = e^^ *(a) = - Te^-da. 

o Jo 

In terms of these matrices, we have 

Ti =(r - mT)^{mT)^{T - mT)G, 
T2 =mr*(mr)G. 



(2.70) 



(2.71) 



The definition (2.70) is also useful from a computational point of view. If 
we recall the series definition of the matrix exponential, 



fc=0 



then we get 



But now we note that the series for $(a) can be written as 
If we let A; = > + 1 in the sum, then, as in (2.58), we have 



= i+T 

= / + a*(o)F. 



aF 



(2.72) 



(2.73) 
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(2.74) 



single sum we can compute * and T ^^'^ *^ 

If we return to the case ^ = 0. m ^ 0. the discrete state equations are 

xik + 1) = *x(*) + Fmik) + r2u{k + 1), 

define a new state. = x(.) - r.„(ib). Then tSl^quatit^^^^^ 
^{k + 1) = x(A: + 1) - r2«(* + 1) 

= *x(*) + FMk) + r2u(k + 1) - r2u(ife + 1), 

i(k + 1) = + TatiC*)] + ritt(*) 

= «e(*) + (*r2 + ri)u(jfe) 
= 4e(fc) + r«(*). 

The output equation is 

yik) = Hx(A:) 

= H{e(A:) + r2u(A;)] 
= He(A;) + Hr2«(ife) 

= U^ik) + J^uik). (2.75) 

Thus for £ = 0, the state equations are given bv (2 VU O ta\ « j 
Note especiaUy that if m = 0 then T, - n «nJ ^' ^ f ^' 
the previous ^odel with no delly ' ~ '^"^ ^"^^'^"^ '^"^ '° 

Our next case is ^ = 1. Rrom (2.67), the equations are given by 

x(* + 1) = #x(A;) + Tiu{k - 1) + T^uik). 
dtaendon of the state, <uk1 tK iiiuol aie ^™ 
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0 be computed,The 
) be and from this 

ate equations are 

■), 

these equations in 
- 1). To do this, we 
tions are 



(2.74) 



(2.75) 

;2.74), and (2.75). 
luations reduce to 

! given by 



nd side, which we 
(hus an increased 



!«(*), 







2-1 • 


p 



2-' 






^ 




V 

/delays 



Figure 2.20 Block diagram of system with delay of more than one period. Double 
line indicates vector valued variables. 

For our final case, we consider ^ > 1. In this case, the equations are 

x(jb + 1) = *x(ifc) + Tiu{k + T2u{k 

and we must eliminate the past controls up to u{k). To do this we introduce 
i new variables such that 

Xn+l(*) = ti(fc - i), Xn+2{k) Xn^i{k) = u{k - 1). 

The structure of the equations is 



■x(* + l) 




aJn+i(fc + l) 




Xn+iik+1) 




.xn+tik + l). 





# Ti r2 0 

0 0 10 

0 0 0 1 

0 0 0: 



»(*) = [H 0] 



X 





\m 1 




"0" 








0 




X„+2{k) 


+ 


0 




.Xn+t(k) . 




.1. 



(2.77) 



(2.76) 



This final situation is easily visualized in terms of a block diagram, as 
shown in Fig. 2.20. 

2*4.5 Numerical Considerations in Computing # and T 

The numerical considerations of these computations are centered in the ap- 
proximation to the infinite sum for ^ given by (2.72) or, for o = T, by (2.58). 
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The problem is that if FT is large, then fPn''/Wlh«™™ - 
before it beconie, small, and £fU "cSable i • 
computer mm>b.r repre^tation, wai^:^^^^ ^h""^ 

a^The basic i.leacom.C(X?C^Tj1.!:^ 



"^iy4 " X::rc~:^?^r,^^ the 

we find a * such that r/2* is f«« i ^ ^ ^n, untU 

;^^^^ue..Wepro^:.l:;r:-- 

we estimate thesise „,R ^ ttesi. .,rher. t^rTI^i: 

fta(Fr)''/(jv+i)i2»» 

Is to select * such that ^ "^J' The rale 

2* > II FT 11 = max 2; 1 1,^. I 

»=1 

Taking the log of both sides, we find 

*>Iog2||Fr|j, 



fi-om which we select 



* = inax(riog2||Fr(|,0). (279) 
imeger and zero is taken because it is possible th«t 11 1:* 



max- 

possible that ||Fr|| is 




s extremely large 
is realized most 
i the value of the 
used by Kahnan 
er and Van Loan 
'and<i-<o = r, 



(2.78) 

I and square the 
and so on, until 
t for deciding on 
J^, which can be 



ample period is 
•m su^ests that 
xl in ^, namely, 



divided by 2* is 
iverge. The rule 



(2.79) 

m X. The max- 
ahat IIFTII is 
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1. Select Fandr. 

2. Conament: Compute || FT ||. 

4. ib 4— smallest nonnegative integer 
greater than log2 V. 

5. Conunent: compute *(r/2*). 

6. Ti ^ r/2* 

7. I ^ Identity 
8. 

9. j ^ 11 

10. If i = 1, 



11. * 



r FTi ^ 



12. i^i-1 

13. Go to step 10. 

14. Conmient: Now double 
9 ib times. 

15. If = 0, stop. 
FT 

16. * " 



17. fc 4- Jfc - 1 
go to step 14. 

18. Go to step 15. 



Figure 2.21 Logic for a program to compute 9 using automatic time scaling. 

already so small that its log is negative, in which case we want to select 
jb = 0. 

Having selected lb, we now have the problem of computing *(r) from 
*(T/2*). Our original concept was based on the series for which satisfied 
(2.78). To obtain the suitable formula for we use the relation between * 
and 9 given by (2.58) as follows to obtam the doubling" formula for 9: 

*(2r) = *(r)*(r), 
I + 2rF*(2r) = [I + rF*(r)][i + rF*(r)] 

= I + 2rF*(r) + T^F^^\T)i 

therefore 

2rF*(2r) = 2rF*(r) + T^F^^^iT). 

This is equivalent to 

*(2r) = (i+ ^*(r)) *(r), 

which is the form to be used. The program logic for computing ^ is 
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Shown in Fig. 2 21.1^ This algorithm does not include the delay discussed 
Fig 2''2^S ^""^ ™Pl«™«°* the logic showH 



2.5 SIGNAL ANALYSIS AND 
DYNAMIC RESPONSE 

In Section 2.3 we demonstrated that if two variables are related by a linear 
constant difference equation, then the ratio of the z-transform of the oZS 

ratio 18 called the transfer function. A method for study of hnear constant 
discrete systems is thereby indicated, consisting of the fiuow^^ter 

1. Compute the transfer function of the system H{z). 

2. Compute the transform of the input signal, E(z). 

3. Form the product, E(z)H{z), which is the transform of the output signal, 

4. Invert the transform to obtain u(kT). 

if the !i^„f ^ ^ difference-equation form, and 

If the mput signal is elementary, then the first three steps of this process 

S;Zd 'V k"* " --Potation. The final step, ho^r. ^^T^l 
done by hand and, because we wiU later be preoccupied with design of trans- 

mtmtion for the kmd of response to be expected from a transform without 
actuaUy invertmg it Our approach to this problem is to presentT^^i^^ 
'm ttT T ''^^^ "1*^ *° their repreTtS 

we will hrlT "f'''^^' '^^^^ Siven an mJcnown transfomi" 

tFhese tomwx. solute, t^imt m^ 

k^^l^ L^m^ ^^'TV'^ *^ whether the im- 

^Zlutilr^o Z ''ir' '° '^^r °f time-response 

t W ^! ^ ^^'^ P™*^"^ attaching a connection between the 

W e^Zt^r^^^^- - ~ ^-^orrns of^I 



^^See X-C2D in Table E.l 
^^'See DELAY in Table E.l. 



2.5 ! 



2.5.1 The Unit 
We have already s 



therefore we have 



This result is muc 
of the unit impuk 
The quantity 
to systems: To cht 
is the imit pulse r 

2.5.2 TheUnil 
Consider the unit 



In this case, the z 



^^We have shifted i 
use subscripts to ic 
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2.5.1 The Unit Pulse 

We have already seen that the unit pulse is defined hy^^ 



ei{k) = 1 
= 0 



(fc = 0) 



therefore we have 



(2.80) 



This result is much like the continuous case, wherein the Laplace transform 
of the unit impulse is the constant 1.0. 

The quantity Ei{z) gives us an instantaneoiis method to relate sign^ils 
to systems: To characterize the system H{z), consider the signal which 
is the unit pulse response; then U{z) == H{z). 

2.5.2 The Unit Step 

Consider the tmit step function defined by 

e2(fc) = l (fc>0) 
= 0 (Jk < 0) 



^l(fc). 



In this case, the z-transform is 



fc=-oo k=0 
1 



z 



z-1 



(|z-M < 1) 

i\z\ > 1). 



(2.81) 



''We have shifted notation here to use e{k) rather than e* for the kth sample. We 
use subscripts to identify different signals. 
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(a) 





Here the transform is characterized by a zero at ^ - 0 anH » « i * 
The significance of the convergence bl^g ^^t'lted Vtlml = 
plored later when we rnnQiHor Tk^ • ^ 1^1 > 1 will be ex- 

transform of t^uIftXt 1 A r'""' '1^^^^ Laplace 

j\ ~ convergence outside the miit circle \z\ - i ™in 

corr^pond to a constant for positive time and zero ^ neS t me 

we sStX^ ^igV 2 Tr^Z:^:^ domailC^r^tplane, 
^(.) marked xffetifirarl^'^'^^^^^ 

time plot of e2[k). ^"P^^^' ^ sketch the 



2.5.3 Exponential 

The one-sided exponential in tune is 



e3(*) = r* (jfc > 0) 

= 0 (*;<0), 



which is the same as r*l(*), using the symbol 1(*) for the 



(2.82) 



unit step function. 
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(b) 



Figure 2.23 (a) Pole and zero of Ez(z) in the z-plane. (b) Plot of esik). 



Now we get 



^ 1 
^ 1 - rz-i 
z 



z — r 



{\rz-'\ < 1) 
(N > |r|). 



(2.83) 



The pole ofEz{z) is at z = r. FVom (2.82) we know that €3 (A;) grows without 
boimd if |r| > 1. Prom (2.83) we conclude that a z-transform that converges 
for large z and has a real pole outside the circle \z\ = 1 corresponds to a 
growing signal. If such a signal were the imit-pulse response of our system, 
such as our digital control program, we would say the program was unstable 
as we saw in (2.37). We plot in Fig. 2.23 the z-plane and the corresponding 
time history of Ez{z) as e3(fc) for the stable value, r = 0.6. 

2.5.4 General Sinusoid 

Our next example considers the modulated sinusoid e4(A;) = [r^cos A:^]1(A;), 
where we assume r > 0. Actually, we can decompose 64 (fc) into the sum of 
two complex exponentials as 



e4(*)=r* ^ 
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and because the z-tramform is linear,^^ we need only compute the transform 
ofeach single complex exponential and add the results later. We thus tZ 

e,{k) = rV^^UA) (2.84) 

and compute 



oo 

-A: 



k=0 

oo 



3=0 



1 



1 - rei^z-i 
z 

= z _ reif (1^1 > (2.85) 

The signal e^{k) grows without bound as k gets large if and only if r > 1 
and a system with this pulse response is BIBO stable if and only if Irl < 1 
The boundary of stability is the unit circle. To complete the argmne Jt Lven 
above for r*cosA:^l(A), we see immediately that 

found by replacing e by -9 in (2.85), ^ 

Z{rVMl(,)} = __^ (N>r). (2.86) 

and thus that 

2\z-reJ»^ z-re-i»j 
_ z(z-rcosd ) 

z^-Mcos&)z + r^ (1^1 >r). f2;87) 

i'n nr2^Tf^;l"^^^^^^^ o^^^C*) - shown 

We note in passing that if ^ = 0, then 64 reduces to eg and, with r = 1 
to e,, so that three of our signals are special case, of e^. By eipToiui tt 

^wTot shown this formaUy. The demonstration, using the definition of 
linearity given above, is simple and is given in Section 2.7. 
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1.0 



(b) 

Figure 2.24 (a) Poles and zeros of E^{z) for 0 = 45**, r = 0.7 in the 2-plane. (b) 
Plot of e4(fc). 



features of E4{z), we can draw a number of conclusions about the relation 
between pole locations in the z-plane and the time-domain signals to which 
the poles correspond. We collect these for later reference. 

1. The settling time of a transient, defined as the time required for the signal 
to decay to one percent of its maximum value, is set mainly by the value 
of the radius, r, of the poles. 

a) r > 1 corresponds to a growing signal that will not decay at all. 

b) r = 1 corresponds to a signal with constant amplitude (which is not 
BIBO stable as a pulse response). 

c) For r < 1, the closer r is to 0 the shorter the settling time. The 
corresponding system is BIBO stable. We can compute the settling 
time in samples, iV, in terms of the pole radius, r. 



pole radius, 


response duration. 


r 


N 


0.9 


43 


0.8 


21 


0.6 


9 


0.4 


5 



d) A pole at r = 0 corresponds to a transient of finite duration. 

2* The number of samples per oscillation of a sinusoidal signal is deter- 
mined by 9. If we require cosOk = cos(^(A: + iV)), we find that a period 
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7* = sampling period 
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j-l samples/cycle. 



unit circle wiS sS n Ji^ts ' ^ ^''^ ^^^^^^^ ^ sketch of the 

Pies per cycle :::^-ZZ:Z%T^^ ThHl'"," 

Hertz is 1/T, and the signal Wnt is / - ^ /^^^ ^T'""'^ 
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2*5.5 Correspondence with Continuous Signals 

jVom the calculation of these few ^-transforms, we have established that the 
duration of a time signal is related to the radius of the pole locations and 
the number of samples per cycle is related to the angle, 6, Another set of 
very useful relationships can be established by considering the signals to be 
samples from a continuous signal, e(*), with Laplace transform E{3), With 
this device we can exploit our knowledge of ^-plane features by transferring 
them to equivalent z-plane properties. For the specific numbers represented 
in the illustration of 64, we take the continuous signal 

y{t) ^ e-^ cosbtl{t) (2.88) 

with 

qT = 0.3567, 
bT = 7r/4. 

And, taking samples one second apart {T = 1), we have 

y(]fer) = (e-^-^«Ycos^l(ib) 

= (0.7)*cos^l(A:) 
= C4(fc). 

The poles of the Laplace transform of y{t) (m the s-plane) are at 

= -a + jh, -a - j6. 
FVom (2.87), the z-transform of E^{z) has poles at 

but b^use y{kT) equals e^^U), ifMlSwsthit 

r = e"*^^, e = 6r, 

If E{z) is a ratio of polynomials in z, which will be the case if e{k) is 
generated by a linear difference equation with constant coefficients, then by 
partial fraction expansion, E{z) can be expressed as a sum of elementary 
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(a) (b) 
Fjure 2,26 Corresponding lines in the 5-plane and the z-plane according to z = 

terms like E4 and E^.^^ In all such cases, the discrete signal can be generated 
by samples from continuous signals where the relation between the 5-plane 
poles and the corresponding z-plane poles is given by 

^ = (2.89) 

T*"^' "^^^"^ ^^^^ ^ P^^^ » place in the 5-plane 

then (2.89) shows us where to look in the z-plane to find a representation of 
discrete samples havmg the same time features. It is useful to sketch several 
major features from the ^-plane to the z-plane according to (2.89) to help 
fix these ideas. Such a sketch is shown in Fig. 2.26. 

Each feature should be traced in the mind to obtain a good grasp of the 
relation. These features are given in Table 2.2. We note in passing that the 
map 2; = e of (2.89) is many-to-one. There are many values of s for each 
valueofz. Infact, if 



then e^i^ = The (great) significance of this fact wiU be explored in 
Chapter 3. 



^^Unless a pole of E(z) is repeated. We have yet to compute the discrete version 
01 a signal corresponding to a higher-order pole. The result is readily shown to be 
a polynomial m ib multiplymg r*e^*^^ 
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ftble 2.2 Descriptions of corresponding lines in 3-plane and r-plane. 



d-plane 



Symbol 



z-plane 



[ Real frequency axis 

5 = <7>0 

5 = <y <0 



= -a + jb 

Constant damping ratio 
if C is fixed and Un 
varies 

s^±j{ir/T) + <T,a<0 



XXX 

□□□ 

OOO 

AAA 



nz| = i 

\ Unit cii 



circle 
z = r > 1 
2 = r, 0<r<l 

z = re-'^ where r = exp(— C'^^n^) 



Logarithmic spiral 



z = -r 



ce signal can be generated 
mtion between the a-plane 
by 

(2.89) 
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»mpute the discrete version j 
jsult is readily shown to be | 



2.5.6 Step Response 

Our eventual purpose, of course, is to design digital controls, and our interest 
in the relation between z-plane poles and zeros and time-domain response 
comes from our need to know how a proposed design will respond in a 
given dynamic situation. The generic dynamic test for controls is the step 
response, and we will conclude this discussion of discrete system dynamic 
response with an examination of the relationships between the pole-zero 
patterns of elementary systems and the corresponding step responses for a 
discrete transfer function from u to y of a hypothetical plant. Our attention 
will be restricted to the step responses of the discrete system shown in 
Fig. 2.27 for a selected set of values of the parameters. 

Note that if zi = pi, the members of the one pole-zero pair cancel 



U{2) 



{2-Z,)(2-Z^) 



nz) 

— o 



Figure 2.27 Definition of the parameters of the system whose step responses are 
to be catalogued. 
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2.! 



I- 




Figure 2.28 Pole-zero pattern of Y{z) for the system of Fig. 2.27, with zi = 
PuU{z) = z/{z - l),oi and a2 selected for 0 = 18*", and C = 0.5. 

the system response, to the input with transform U{z) = 1 (a xmit 
pulse) is 



Y{z) = 



z — rcosO 



z2_2rcosfiz + r2' 



(2.90) 



This transform, when compared with the transform E4{z) given in (2.87), 
is seen to be 

Y{z) = z-'E4{z), 

and we conclude that under these circiunstances the system pulse response 
is a delayed version of 64 (fc), a typical second-order system pulse response. 

For our first study we consider the effect of zero location. We let zi = pi 
and explore the effect of the (remaining) zero location, ^2? on the step- 
response overshoot for three sets of values of ai and 02. We select ai and 
02 so that the poles of the system correspond to a response with damping 
ratio C = 0.5 and consider values of 0 of 18, 45, and 72 degrees. In every 
^^^^^U^a^^^e^ain^^^^^cl^^a|^^e^ 
ec^ilals tn^t^^^^T^^^twIionri^^^^'^^ 

g = 18^. IShewGiu^ lomi^^fe^^ In^i^i^^!^ tp 

the two poles and one zero of H{z)y we show the pole at z = 1 and the zero 
at z = 0, which come from the transform of the input step, U{z)y given by 
z/{z - 1). 

The major effect of the zero Z2 on the step response y{k) is to change 
the percent overshoot, as can be seen from the four step responses for this 
case plotted in Fig. 2.29. To summarize all these data, we plot the percent 
overshoot versus zero location in Fig. 2.30 for C = 0.5 and in Fig. 2.31 



Figure 2.29 Pl< 
pattern of Fig. 2 
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little influence 
it comes near 



1 
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I 
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Figure 2.30 
6 = 18°, 45% 



••27, with zx ^ i 

) = 1 (a unit 

(2.90) 
Mvenin (2.87), 
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o Zj = 0,8 



X z, = 
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0:6 



X • • 



o 

□ 

X o 



y I » < ' 



□ 8 



it- 
Figure 2.29 Plot of step responses for a discrete plant described by the pole-zero 
pattern of Fig. 2.28 for various values of Z2- 

for C = 0.707. The major feature of these plots is that the zero has very 
little influence when on the negative axis, but its influence is dramatic as 
it comes near +1. Also included on the plots of Fig. 2.30 are overshoot 



1000 r- 




1.0 



1.5 



0 0.5 
Zero location Zj 

Figure 2.30 Effects of an extra zero on a discrete second-order system, C = 0.5: 
^ = 18% 45^, and 72**. 



II 
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3 - 



-1.0 0.0 1,0 

Zero location 

Figure 2.31 Effects of extra zero on second-order system when C = 0.707; 0 = 18** 
45°, 72°. Percent overshoot versus zero location. 



figures for a zero in the unstable region on the positive real axis. These 
responses go in the negative direction at first, and for the zero very near +1, 
the negative peak is larger than one!^^ 

Our second class of step responses corresponds to a study of the influence 
of a third pole on a basically second-order response. For this case we again 
consider the system of Fig. 2.27, but this time we fix zi = Z2 = -1 and 
let pi vary from near ~1 to near -hi. In this case, the major influence of 
the moving singularity is on the rise time of the step response. We plot this 
fSrp^^^i Ar ^^N& ^ ^w? w ^ P i tt T i^ "^ 9; ^ ?^ 2,32, In the figure 

whjcjyaio^jpOj^ • - ' 

rise time to get very nduch longer as U 
and comes to dominate the response. 



Such systems are called nonminimum phase by Bode because the phase shift they 
impart to ^ sinusoidal input is greater than the phase of a system whose magnitude 
response is the same but that has a zero in the stable rather than the unstable 
region. 




lOOi- 
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-1.0 0.0 

Pole location 

Figure 2.32 Effects of extra pole on system rise time- Two zeros at -1, one zero 
at oo; C = 0.5; 6 = 18^ 45% 72**. 



Our conclusions from these plots are that the addition of a pole or zero 
to a given system has only a small effect if the added singularities are in the 
range from 0 to -1. However, a zero moving toward z = +1 greatly increases 
the system overshoot. A pole placed toward z = +1 causes the response to 
slow down and thus prhnarily affects the rise time, which is progressively 
increased. 



2,6 FREQUENCY RESPONSE 

A very important concept in linear systems analysis is the frequency re- 
sponse. If a sinusoid at frequency u)o is applied to a stable, linear, constant, 
continuous system, the response is a transient plus a sinusoidal steady state 
at the same frequency, Uo, as the input If the transfer function is written 
in gain-phase form as H{j(jj) = >l(a;)eJ'^*^), then the steady-state response^ 
to a luiit-amplitude sinusoidal signal has amplitude A{(Jo) and phase V'(^o) 
relative to the input signal. 

We can say almost exactly the same respectmg the frequency response 
of a stable, linear, constant, discrete system. If .the system has a transfer 
function H{z), we define its magnitude and phase for z taking on values 
around the unit cfrcle by /r(e^*^^) = i4(a;r)eJ**<^^>. If a unit-amplitude 
sinusoid is applied, then in the steady state, the response samples will be on 
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a smusoid of the same frequency with amplitude Aiu,,T) and phase ^(u,„T) 
It worthwhile going through the calculations to fix ideas on this point 
Worn (2.16), the discrete response transform is 



(2.91) 
■ (*>oT, we have 

(2.92) 



(2.93) 



U(z) = Hiz)E{z). 
If e(k) = cos(u;oTk)l(k), then, from (2.87) with r = 1 and ^ 

Eiz) = i I £ + — i_\ 

2 U - ej"oT ^ g_ e-j'-'oTf ' 

If we substitute (2.92) into (2.91), we obtain 

C^(z) = l|^^ + __f^(fL\ 

The steady state of u{kT) corresponds to the terms in the expansion of 
(2.93) associated with the two poles on the unit circle If we exZd7/?.\/ 

£/,,(z) = l£(£!rf)f , l H(e-i-oT^z 
2 z- ei'^oT ^ 2 z- e-j"oT ' 

If H(e^oT^ = ^(c.o2')e'>(«<'^, then we have 

.4 e^'^z 



and the inverse transform of Ugs{z) is 



+ 4 



2 z-e-i'^o^' 



(2.94) 



2 2 
= >lcos(a;„I'* + V'), (2.95) 

rh» "^1^ of particular frequency responses untU later 

chapters (see, for example, Figs. 4.3. 4.8, 5.16, and 5.23). Ho^r It sho,S 

thro k1 T "'^^^"^^ * '"^"^'^ of frequency "^ li*ptS 
quency a,, + ^^/r for mteger £ which also pass through these points. This 
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is the phenomenon of aliasing, to which we will return in Chapter 3. Here, 
we repeat, we define the discrete frequency response of a transfer function 
H{z) to sinusoids of frequency ojo as H{€?^^'^). 



The Discrete Fourier IVansform (DFT), The analysis developed 
above based on the z-transform is adequate for considering the theoreti- 
cal frequency response of a linear, constant system or the corresponding 
difference equation, but it is not the best for the analysis of real-time signals 
as they occur in the laboratory or in other experimental situations. For the 
analysis of real data, we need a transform defined over a finite data record, 
which can be computed quickly and accurately. The required formula is 
that of the Discrete Fourier Transform, the DFT, and its numerical cousin, 
the Fast Fourier IVansform, the FFT. Implementation of a version of the 
FFT algorithm is contained in all signal-processing software and in most 
computer-aided control-design software. 

To understand the DFT, it is useful to consider two properties of a 
signal and its Fourier transform that are complements of each other: the 
property of being periodic and the property of being discrete. In ordinary 
Fourier analysis, we have a signal that is neither periodic nor discrete and 
its Fourier transform is also neither discrete nor periodic. If, however, the 
time function f{t) is periodic with period To, then the appropriate form 
of the transform is the Fourier series, and the transform is defined only 
for the discrete frequencies u) == 27rn/ro. In other words, if the function 
in time is periodic, the function in frequency is discrete. The case where 
the properties are reversed is the z-transform we have just been studying. 
In this case, the time functions are discrete, being sampled, and the z- 
transform is periodic in u)\ for if z = eP'^'^ , corresponding to real frequencies, 
then replacing c*; = a; + 27r k/T leaves z unchanged. We can summarize these 
results with the following table: 






Time 


Frequency 


Fourier series 


periodic 


discrete 


z-transform 


discrete 


periodic 



Suppose we now have a time function that is both periodic and discrete. 
Based on what we have seen, we would expect the transform of this function 



80 CHAPTER 2 SYSTEMS ANALYSIS 



also to be both periodic and discrete. And this is the case, which leads 
us to the finite discrete Fourier transform and its finite inverse. Let the 
time function in question be f{kT) = f{kT + NT). Because the function is 
periodic, the transform can be defined as the finite sum 

^^^^ ^ ifc=0 

This is the same as the z-transform over one period evaluated at the dis- 
Crete frequencies of a Fourier series a; = 2nn/NT. It is standard practice 
to suppress aU the arguments except the indices of time and fi-equency and 
write 



fc=0 



(2.96) 



To complete the DFT, we need the inverse transform, which, by analogy 
with the standard Fourier transform, we guess to be the sum 

N-l 
n=0 

If we substitute (2.96) with summing mdex ^ into this, we find 

JV-l (N-l 



n=0 W=0 

Intercfaan^ng the order of the summations gives 



t=Q I n=0 J 



The sum in the braces is a finite geometric series, which we can evaluate as 
follows: 

3r e^-2»[n(ik-/)]/Ar ^ 1 - eJ^'^Cfc-O 
^ 1 _ ei'^Ak-D/N 



■I 



N k-£=0 

0 A:-^=l,2,...,JNr-L 



The sum is pi 
sum we have i 
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Because there 
compute the ] 
tiply and add 
it would take 
several author 
take advantagt 
of Fn can be * 
a power of 2. 
large value. Tl 
IVansform, or 
To use the 
have a system 
quency uji = 
input to the s 
time, the outp 
is 



i 



2.6 FREQUENCY RESPONSE 81 



iich leads 
Let the 
action is 



• the dis- 
practice 
2ncy and 



(2.96) 
analogy 



The sum is periodic with period N. With this evaluation, we see that the 
sum we have been considering is Nfkj and thus we have the inverse sum, 



n=0 

Equations (2.96) and (2.97) comprise the DFT: 



(2.97) 



N-i 



k=0 



n=0 

Because there are N terms in the sum in (2.96), it would appear that to 
compute the DFT for one frequency it will take on the order of N mul- 
tiply and add operations, and to compute the DFT for all N frequencira, 
it would take on the order of multiply and add operations. However, 
several authors, especially Cooley and Tukey(1965), have showed how to 
take advantage of the circular natiure of the exponential so that all N values 
of F„ can be computed with on the order of N log(JV) operations if N is 
a power of 2. For N = 1024, this is a saving of a factor of 100, a very 
large value. Their algorithm and related schemes are called the Fast Fourier 
Ti-ansform, or FFT. 

To use the DFT/FFT in evaluating frequency response, we suppose we 
have a system described by (2.91) and that the input is a sinusoid at fre- 
quency ut = 2ire/NT so that e(kT) = Asm{2nikT/NT). We apply this 
input to the system and wait untU all transients have died away. At this 
time, the output is given by u(kT) = B 8m{2irtk/N + V). The DFT of e(fc) 
is 



£„=£Xsin(H^)e-«'' 

1 n \ / 



fc=0 

ik=0 

0, e^h, 

= < NA 
2i' 



£ = n. 
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The DFT of the output is 

N-l 



k=Q ^ ^ * 

fo. 



^eP^ i = n 

Dividing these results, we see that with sinusoidal input and output, the 
frequency response at the fi^uency u = {2it£)/NT is given by 

El 

where Ui = FFT{uk) and Et = FFT{ek), each evaluated at n = ^. We will 
discuss in Chapter 8 the general problem of estimation of the total frequency 
response from experimental data usmg the DFT/FFT as well as other tools. 



2.7 PROPERTIES OF THE ^-TRANSFORM 

We have used the z-transform to show that Imear, constant, discrete systems 
can be described by a transfer function that is the z-transform of the system's 
unit-pulse response, and we have studied the relationship between the pole- 
zero patterns of transfer functions in the z-plane and the corresponding time 
responses. We began a table of z-transforms, and a more extensive table is 
given m Appendbc B. In Section 2.7.1 we turn to consideration of some of the 
properties o f the z -transform that are essential to the effective and correct 

tramfer fi mction is Min. 

2.7.1 z-IVansform Properties 

In order to make maxunum use of a table of z-transforms, one must be able 
to use a few simple properties of the z-transform which follow directly from 
the definition. Some of these, such as linearity, we have afready used without 
making a formal statement of it, and others, such as the transform of the 
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convolution, we have previously derived. For reference, we will demonstrate 
a few properties here and collect them into Appendix B for future reference. 
In all the properties listed below, we assume that Fi{z) = Z{fi{kT)}. 

1. Linearity: A function f{x) is linear if f{axi + PX2) = af{xi) + Pf{x2). 
Applying this result to the definition of the z-transform, we find imme- 
diately that 

k=—oo 

^aZ{h{k)}+0Z{f2{k)} 
= aFi{z) + pF2{z). 

Thus the z-transform is a linear function. It is the linearity of the trans- 
form that makes the partial-fraction technique work. 

2* Convolution of Time Sequences: 



■2 1^ E /i(0/2(* - 0 j = Fi(z)F2{z). 

We have already developed this result in connection with (2.30). It is 
this result with linearity that makes the transform so useful in linear- 
constant-system analysis because the analysis of a combination of such 
dynamic systentis can be done by linear algebra on the transfer functions. 

3. Time Shift: 

Z{f{k + n)} = z+"F(z). (2.98) 
We demonstrate this result also by direct calculation: 

Z{f{k + n)}= f; f{k + n)z-^. . 
ik=-oo 

If we let + n = j, then 

2{/(fc+n)}= f; /o>-o-«) 

j=-oo 

= QED 
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This property is the essential tool in solving linear constant-coefficient 
difference equations by transforms. We should note here that the trans- 
form of the time shift is not the same for the one-sided transform because 
a shift can introduce terms with negative argument which are not in- 
cluded in the one-sided transform and must be treated separately. This 
effect causes initial conditions for the difference equation to be intro- 
duced when solution is done with the one-sided transform. See Problem 
2.13. 



4. Scaling in the z-plane: 



(2.99) 



By direct substitution, we obtain 



Final' Value 



The conditio 
strictly insid 
in {z " l)F{ 
nitude of z \ 
k. Therefore 
the possible 
The size of 1 
partial-fracti 



r 



fc=-oo 
= E /(fc)(rz)-* 

A:=— OO 

= F{rz). QED 

As an illustration of this property, we consider the z-transform of the 
imit step, l(fc), which we have computed, before: 



-k 



2{i(*)} = E^ 

k=0 



By property 4 we have immediately that 



z-1 



Z{r-H{k)} = 



rz 



As a more general example, if we have a polynomial a{z) = z^ + aiz + 
02 with roots re^-^^, then the scaled polynomial a^z^ + aiaz + 02 has 
roots {r/a)e^^^. This is an example of radial projection whereby the 
roots of a polynomial can be projected radially simply by changing the 
coefficients of the polynomial. The technique is sometimes used in pole- 
placement designs as described in Chapter 6, and sometimes used in 
adaptive control as described in Chapter 11. 
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5, Final'Value Theorem: If F{z) converges for \z\ > 1 and all poles of 
{z — l)F(z) are inside the unit circle, then 



6. 



Um f{k) = ]im{z ~ l)F(z). 



(2.100) 



The conditions on F{z) assure that the only possible pole of F{z) not 
strictly inside the unit circle is a simple pole at z = 1, which is removed 
hi (z - l)F{z). Furthermore, the fact that F{z) converges as the mag- 
nitude of z gets arbitrarily large ensures that f{k) is zero for negative 
k. Therefore, all components of f{k) tend to zero as k gets large, with 
the possible exception of the constant term due to the pole at z = 1. 
The size of this constant is given by the coefficient of l/(z - 1) m the 
partial-fraction expansion of F(z), namely, 

C = lim(z - l)F{z). 

However, because all other terms in f{k) tend to zero, the constant C 
is the final value of f{k), and (2.100) results. QED 

As an illustration of this property, we consider the signal whose trans- 
form is given by 



U{z) = 



z Tz+1 



\z\ > 1. 



Z-0.5 2Z-1' 
Because U{z) satisfies the conditions of (2.100), we have 

lim u(k) = lim(z — 1) — -;rT?~~~T 
jb-»oo ^ ' z-0.5 2 z-1 



= lim 



« - 0.5 2 



(^ + 1) 



1 T 



1 - 01^2 

= 2r. 



This result can be checked against the closed form for u{k) given by 
(2.106) below. 

Inversion: As with the Laplace transform, the z-transform is actually 
one of a pair of transforms that connect functions of time to functions of 
the complex variable z. The z-transform computes a function of z from 
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a sequence in k. (We identify the sequence number A; with time in our 
analysis of dynamic systems, but there is nothing in the transform per 
se that requires this.) The inverse z- transform is a means to compute a 
sequence in Jfc from a given function of z. We first examine two elementary 
schemes for inversion of a given F{z) which can be used if we know 
beforehand that F{z) is rational in z and converges as z approaches 
infinity. For a sequence /(fc), the z-transform has been defined as 

n^)^ E /(*)^"*' ro<\z\<Ro. (2.101) 

If any value of f{k) for negative k is nonzero, then there will be a term 
in (2.101) with a positive power of z. This term will be unbounded if 
the magnitude of z is imbounded; and thus if F{z) converges as |z| 
approaches infinity, we know that f{k) is zero for A: < 0. In this case, 
(2.101) is one-sided, and we can write 

i^W = E ^0<I4 (2.102) 

k=0 

The right-hand side of (2.102) is a series expansion of F(z) about infinity 
or about z~^ =0. Such an expansion is especially easy if F{z) is the 
ratio of two polynomials in z~^. We need only divide the nimierator by 
the denominator in the correct way, and when the division is done, the 
coefficient of z"* is automatically the sequence value f{k). An example 
we have worked out before will illustrate the process. Suppose we take 
our system to be the trapezoid-rule integration with transfer function 
given by (2.14): 

We will take the input to be the geometric series represented by e3(A;) 
with r = 0.5. Then we have 
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U{z)^E3{z)H{z) 

_ z T z+l 
z-0.5 2"z-l' 



\z\ > L 



(2.103) 



Equation (2.103) represents the transform of the system output, u{k). Keep- 
ing out the factor of r/2, we write U{z) ssa ratio of polynomials in z'^. 



U{z) = 



1 + z 



-1 



2 l-L5z-i+0.5z-2' 



(2.104) 



and divide as follows: 



^[1 + 2.5z'^ + 3.25z~^ + 3.625z-^ + • • 

St 

1 - 1.5z-^ +0.5z-^)l + z"^ 

1- 1.5z~^ + 0.5z~^ 
2.5z-^ - 0.5z-^ 
2.5z-^ - 3.75z~' + 1.25z~^ 
3.252~2 - 1.25z-^ 
3.25z~^ - 4.875z~^ + 1.6252-" 
3.625z-^ - 1.625Z-'* 
3.625z-^ 

By direct comparison with U{z) = ]Co° «(*')^~*» conclude that 

«o = r/2, 

ui = (r/2)2.5, 

«2 = (T/2)3.25, (2.105) 



' ies represented by ejik) 



Clearly, the use of a computer will greatly aid the speed of this process in all 
but the simplest of cases. Some may prefer to use synthetic division and omit 
copying over all the extraneous z's in the division. The process is identical 
to converting F{z) to the equivalent difference equation and solving for the 
unit-pulse response. 
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The second special method for the inversion of z-transforms is to de- 
compose F(z) by partial-fraction expansion and look up the components of 
the sequence /(*) in a previously prepared table. We consider again (2.103) 
and expand U{z)asa function of as follows: 



Uiz) = 



B 



2 1 - z-i 1 - O.Sz-i 1 - z-i 1 _ 0.5Z-1 
We multiply both sides 1 - z"!, let z"^ = 1, and compute 



T 2 
2 0.5 



Sunilarly, at z~* = 2, we evaluate 



B = 



T l + 2 
2 1-2 



3T 



Looking back now at C2 and es, which constitute our "table" for the moment, 
we can copy down that 

nrp 



(-fG)>(*) 



(2.106) 



Evaluation of (2.106) for * = 0, 1, 2, . . . wiU, naturally, give the same values 
for «(*) as we found in (2.105). We now examine more closely the role of the 
region of convergence of the z-transform and present the inverse-transform 
integral. We begin with another example. The sequence 



/(*) 



1 0, 



A;<0, 
*>0, 
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has the transform 



-1 



A:=— oo 



\z\ < 1. 



This transform is exactly the same as the transform of the unit step 
1(A:), (2.81), except that this transform converges inside the unit circle and 
the transform of the l{k) converges outside the unit circle. Knowledge of 
the region of convergence is obviously essential to the proper inversion of 
the transform to obtain the time sequence. The inverse z-transform is the 
closed, complex integral^^ 



/(*) 



.dz 



(2.107) 



where the contour is a circle in the region of convergence of F{z). To demon- 
strate the correctness of the integral and to use it to compute inverses it is 
useful to apply Cauchy's residue calculus [see Churchill and Brown (1984)]. 
Cauchy's result is that a closed integral of a function of z which is ana- 
lytic on and inside a closed contour C except at a finite number of isolated 
singularities Zi is given by 



^J^F{z)dz = '£Res{z,), 



(2.108) 



In (2.108), Res(zi) means the residue of F{z) at the singularity at Zj. We 
will be considering only rational functions, and these have only poles as 
singularities. If F{z) has a pole of order n at zi, then {z — zi)^F{z) is 
regular at zi and can be expanded in a Taylor series near zi as 



{z -.zi)^F(z) =A_„ + A^n^iiz - zi) + 
The residue of F{z) at zi is A^i. 



+ i4-i(z-zi) 



n-l 



(2.109) 



^^If it is known that f{k) is causal, that is, f{k) = 0 for A: < 0, then the region 
of convergence is outside the smallest circle that contains all the poles of F{z) for 
rational transforms. It is this property that permits inversion by partial-fraction 
expansion and long division. 
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r 



First we will use Cauchy's formula to verify (2.108). If F{z) is the ^- 
transform of /(fc), then we write 



*= — OO 



We assume that the series for F{z) converges uniformly on the contour of 
integration, so the series can be integrated term by term. Thus we have 

dz 



The argument of the integral has no pole inside the contour if ik ~ / > l, 
and it has zero residue at the pole at z = 0 if A: - Z < 0. Only if Jfe = / does 
the integral have a residue, and that is 1. By (2,108), the integral is zero if 
k # / and is 27rj if A; = /. Thus I = /(fc), which demonstrates (2.107). 

To illustrate the use of (2.108) to compute the inverse of a z-transform, 
we will use the function z/{z — 1) and consider first the case of convergence 
for \z\ > 1 and second the case of convergence for \z\ < 1. For the first case, 



fiik) = 



dz 



2irj J\z\=R>i z-1 z 



(2.110) 



where the contour is a circle of radius greater than 1, Suppose Jb < 0. In this 
case, the argument of the integral has two poles inside the. contour: one at 
i = 1 with residue 



and one pole at z = 0 with residue found as in (2. 109) (if A; < 0, then 
removes the pole): 



1 



z-1 



1-z 



= -(1 + z + + • • • + + •••)• 
The residue is thus —1 for all A:, and the sxmi of the residues is zero, and 

/i(fc) = 0, k<0. (2.111) 
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For ik > 0, the argument of the integral in (2.110) has only the pole at z = 1 
with residue 1. Thus 



/i(Jfc) = l, k>0. 



(2.112) 



Equations (2.108) and (2.109) correspond to the unit-step function, as they 
should. We would write the inverse transform symbolically Z~^{,} as, in 
this case, 



(2.113) 



when z/{z — 1) converges for \z\ > 1. 

If, on the other hand, convergence is inside the unit circle, then for A: > 0, 
there are no poles of the integrsmd contained in the contour, and 

f2{k) =0, k>0. 

At k < Oy there is a pole at the origin of z, and as before, the residue is 
equal to —1 there, so 

f2{k) = -1, jb < 0. 
In symbols, corresponding to (2.113), we have 

when z/{z — 1) converges for \z\ < 1. 

Although, as we have just seen, the inverse integral can be used to com- 
pute an expression for a sequence to which a transform corresponds, a more 
effective use of the integral is in more general manipulations. We consider one 
such case that will be of some interest later. First, we consider an expression 
for the transform of a product of two sequences. Suppose we have 

f3ik) = fiik)f2ik), 

and /i and /a are such that the transform of the product exists. An expres- 
sion for Fz{z) in terms of Fi{z) and F2iz) can be developed as follows. By 
definition 



F3{Z)= E /l(*)/2(*)^ 



'1 
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From the inversion integral, (2.107), we can replace f2{k) by an integral: 



*=-oo 



We assume that we can find a region where we can exchange the summation 
with the integration. The contour will be called in this case: 

The sum can now be recognized as Fi{z/Q and, when we substitute this, 

the contour C3 must be in the overlap of the convergence regions of F2{Q 
and Fi{z/Q. Then Fz{z) will converge for the range of values of z for which 
C3 can be found. 

If we let /i = /2 and z = 1 in (2.114), we have the discrete version of 
Parseval's theorem, where convergence is on the unit circle: 

^3(1) = tj = 2iji^.Km (J) f • 

This particular theorem shows how we can compute the sum of squares 
of a time sequence by evaluating a complex integral in the z-domain. The 
result is useful in the design of systems by least squares. 

2.7.2 Another Derivation of the Transfer Function 

Let 2> be a discrete system which maps an input sequence, {e(A:)}, into an 
output sequence {tt(A;)}.^^ Then, expressing this as an operator on e(fc), we 
have 

u{k) = V{e{k)}. 

If P is linear, then 

^>{aei(ik) + I3e2{k)} = aV{ei{k)} + jW>{e2(fc)}- (2.116) 



^^This derivation was suggested by L. A. Zadeh in 1952 at Columbia University. 
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If the system is constant, a shift in e{k) to e{k + j) must result in no other 
effects but a shift in the response, u. We write 



V{e{k + j)} = u{k + j) for all j 



V{e{k)} = u{k). 



(2.117) 



Theorem. If 2? is linear and constant and is given an input for a value 
of z for which the output is finite at time then the output will be of the 
form H{z)z^. 




Proof. In general, if e(fc) = z*, then an arbitrary finite response can be 
written 

u{k)^H{z,k)z^. 

Consider 62(4) = z*^^ = z^z* for some fixed j. Prom (2.116), if we let 
a = z^ it must follow that 

U2 = z^u{k) 
^ziH{z,k)z^ 

= H{zik)z^'^K (2.118) 

Rrom (2.117), we must have 

U2{k) = u{k + j) 

= Jf (z, i + jk)z**^ for all j. (2.1191 

Rrom a comparison of (2.118) and (2.119), it follows that 

H{z, k) = H(z, k + j) for all j; 

that is, H does not depend on the second argument and can be written H{z), 
Thus for the elemental signal e{k) = z*, we have a solution u{k) of the same 
(exponential) shape but modulated by a ratio fi^(z), u{k) = H{z)z^. 
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Can we represent a general signal as a linear sum (integral) of such 
elements? We can, by the inverse integral derived above, as follows: 



where 



(2.120) 



(2.121) 



for signals with r < R tor which (2.121) converges. We call Eiz) the z- 
transform of c(jfc), and the (closed) path of integration is in the annular 
region of convergence of (2.121). If e(A:) = 0, fc < 0, then i? oo, and this 
region is the whole z-plane outside a circle of finite radius. 

The consequences of linearity are that the response to a sum of signals 
is the sum of the responses as ^ven in (2.116). Although (2.120) is the limit 
of a stmi, the result still holds, and we can write 

u(*;) = ^ £(2) [response to z*]— , 
but, by the theorem, the response to z* is H{z)z''. Therefore we can write 
m = -^jE{z)[H{z)A^' 



= ^/./f(z)^(z)z*f. 



(2.122) 



We can define U{z) = H{z)E{z) by comparison with (2.120) and note that 



U{z)^ E u{k)z-^^H{z)E{z). 



(2.123) 



*=-oo 



Thus H{z) is the transfer function, which is the ratio of the transforms of 
e{k) and u{k) as well as the amplitude response to inputs of the form z*. 

This derivation begins with linearity and stationarity and derives the 
z-transform as the natural tool of analysis from the fact that input signals 
in the form produce an output that has the same shape.^^ It is somewhat 
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a) u(A:) = ( 
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more satisfying to derive the necessary transform than to start with the 
transform and see what systems it is good for. Better to start with the 
problem and find a tool than start with a tool and look for a problem. 
Unfortunately, the direct approach requires extensive use of the inversion 
integral and more sophisticated analysis to develop the main result, which 
is (2.123). Chacun k son gout. 



2.8 SUMMARY 

In this chapter we have shown how systems described by linear difference 
equations with constant coefficients can be described by transfer functions 
if the signals are represented by ^-transfornas. The transfer function was 
shown to be the z-transform of the unit-pulse response of the system; and, 
furthermore, the system output was shown to be the convolution of the input 
with the unit-ptdse response. Using this result, we showed a condition for 
the Bounded-Input-Bounded-Output stability of a linear, constant system. 
The test developed by Jury for a polynomial to have all roots inside the unit 
circle was introduced to provide a stability test. We introduced the observer 
and the control canonical forms for transfer functions and gave rules for 
block-diagram reduction of transfer functions. We also introduced the state 
descriptions of these canonical forms and showed how to derive the matrices 
of the dynamic system in state form. We then derived the discrete transform 
for a sampled-data system both by transform and by state-space methods. 
The latter are especially well suited for computer implementation. 

We studied the dynamic response of discrete systems, including espe- 
cially the step response of a second-order system. The effects of the location 
of a zero and of a third pole were plotted, largely for future reference in 
design. 

Several of the properties of the z-transform were demonstrated, and the 
calculation of the inverse of a z-transform was presented by long division, 
by p2urtml-£raction expansion, and by evaluation of the inverse transform 
integral. 



PROBLEMS AND EXERCISES 

2.1 Check the following for stability: 

a) tt(ik) 0.5u(Jb - 1) - 0.3u{* - 2) 

b) u{k) = h6u{k - 1) - u{k - 2) 

c) u{k) = 0.8t4(i - 1) -h OAu{k - 2) 



96 CHAPTER 2 SYSTEMS ANALYSIS 



2.2 a) Derive the difference equation corresponding to the approximation of inte- 

gration found by fitting a parabola to the points efc-2j ^k-u and taking 
the area under this parabola between t = kT — T and * = fcT as the 
approximation to the integral of e{t) over this range, 
b) Find the transfer function of the resulting discrete system and plot the 
poles and zeros in the z-plane. 

2.3 Verify that the transfer function of the system of Fig. 2.8(c) is given by the 
same H{z) as the system of Fig. 2.9(c). 

2.4 a) Compute and plot the unit-pulse response of the system derived in Prob- 

lem 2.2. 

b) Is this system BIBO stable? 

2.5 Consider the difference equation 

ti(A; + 2) = 0.25u(ik). 

a) Assiune a solution ti(ife) = i4<z* and find the characteristic equation in z, 
h) Find the characteristic roots zi and Z2 and decide if the equation solutions 

are stable or unstable. 

c) Assume a general solution of the form 

u(ik) = i4iZi* + i42Z2* 

and find Ai and A2 to match the initial conditions u{0) = 0,u(l) = 1. 

d) Repeat parts (a), (b), and (c) for the equation 

uik + 2) = ~0.25u(fc) 

e) Repeat parts (a), (b), and (c) for the equation 

u(A; + 2)=u(ik + l)-0.5tt(*). 

2.6 Show that the characteristic equation 

, z^-2rcos(e)z + r* 

has the roots 

2.7 a) Use the method of block-diagram reduction, applying Figs 2.5, 2.6, and 

2,7 to compute the transfer function of Fig 2.8(c). 

b) Repeat part (a) for the diagram of Fig. 2.9(c). 
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2.8 Apply Jury's test to determine it fk- r 11 • 

any roots outside the unjf circle "^'"^ characteristic equations have 

a) + 0.25 



b) 
c) 



2^ 0.012 + 0.405 
2' -3.6^2 + 42 -1.6 



^^C^ .ook -up .he ^ ^ 
b) ' 
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the poles and zeros m the z-plane r - 0.5 and plot the location of 

M. U»^. CAD .oo, .. co„p„. ^ ^ ^^^^ ^ 

the results in the 2-plane the zeros and poles of 

ta tbe z-plane. ° ' ™ "» ««• and pote, of the resulte 

-^"-i "y the teu^ 



(« + !)(« + 3)- 



matrices F,G,H, J. ^' «>"espomiing description 

b) in partial fractions and draw the rorr 

diagram with each component p^ttTn!? f^^P^'^dmg parallel block 
state e and give the cor^esS^Vtl^T -^^^^^^ ^"^^ 'be 

c) By finding the transfer ^Zo^'xJut^Tm T^^'J 

fraction form, express x, and ^ * ^^^^ P^* (») Partial 

"-.i.^ « .he r t.f^ 
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d) Verify that the matrices you have found are related by the formulas 

A = T-^FT, 
B = T-^G, 
C = HT, 
D = J. 

2«12 The first-order system (z - a)/(l - a)z has a zero at z = a. 

a) Plot the step response for this system for a = 0.8, 0.9, 1.1, 1.2, 2. 

b) Plot the overshoot of this system on the same coordinates as those ap- 
pearing in Fig. 2.30 for -1 < a < 1. 

c) In what way is the step response of this system unusual for a > 1? 

2.13 The one-sided z-transform is defined as 



-lb 



0 

a) Show that the one-sided transform of f{k + 1) is 

Z{/(*-hl)} = zF(z)-z/(0). 

b) Use the one-sided transform to solve for the transforms of the Fibonacci 
ntunbers by writing (2.4) as tijb+2 = ujb+i + Uk- Let uq = Ui = 1. [ You 
will need to compute the transform of f{k + 2).] 

c) Compute the location of the poles of the transform of the Fibonacci num- 
bers. 

d) Compute the inverse transform of the numbers. 

e) Show that if Uj^ is the ibth Fibonacci number, then the ratio ujb+i/^fc will 
go to (1 + v^)/2, the golden ratio of the Greeks. 

f) Show that if we add a forcing term, e{k), to (2.4) we can generate the 
Fibonacci numbers by a system that can be analyzed by the two-sided 
transfSnn; i.e., let iik = + + and let ejk == Sail^^pc^ = 
1 at ib = 0 and zero elsewhere]. Take the two-sided transform and show 
that the same U{z) results as in part (b). 

2.14 Substitute u = Az^ and e = Bz^ into (2.2) and (2.7) and show that the 
transfer functions, (2.15) and (2.14), can be found in this way. 

2.15 Consider the transfer function 



H{z) = 



(z + l)(z^-1.3z4-0.81) 
(z2 - 1.2z -h 0.5)(z2 - 1.4z + 0.81) ' 
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Draw a cascade realization, using observer canonical forms for second-order blocks 
and in such a way that the coefficients as shown in H{z) above are the parameters 
of the block diagram. 

2.16 a) Write the H{z) of Exercise 2.15 in partial fractions in two terms of second 

order each, and draw a parallel realization, using the observer canonical 
form for each block and showing the coefficients of the partial-fraction 
expansion as the parameters of the realization, 
b) Suppose the two factors in the denominator of H{z) were identical (say 
we change the 1.4 to 1.2 and the 0.81 to 0.5). What would the parallel 
realization be in this case? 

2.17 Show that the observer canonical form of the system equations shown in Fig. 
2.9 can be written in the state-space form as given by (2.26). 

2.18 Draw out each block of Fig. 2.10 in (a) control and (b) observer canonical 
form. Write out the state-description matrices in each case. 

2.19 For a second-order system with damping ratio 0.5 and poles at an angle in 
the z-plane of 9 — 30**, what percent overshoot to a step would you expect if the 
system had a zero at Z2 — 0 6? 

2.20 Consider a signal with the transform (which converges for \z\ > 2) 



(z-l)(z-2)- 



a) 
b) 
c) 

2.21 a) 



b) 



What value is given by the formula (final- value theorem) of (2.100) applied 
to this U{z)l 

Find the final value of u(A:) by taking the inverse transform of {/(z), using 
partial-fraction expansion and the tables. 
Explain why the two results of (a) and (b) differ. 

Find the z-transform and be sure to give this region of convergence for 
the signal 

u(ife) = r+l*l, r<l. 

[Hint: Write u as the sum of two functions, one for fc > 0 and one for 
k < O^^ffind, tfaeindividual transforms, and determine values of z for which 
both terms converge.] 

If a rational function U{z) is known to converge on the unit circle \z\ = 1, 
show how partial-fraction expansion can be used to compute the inverse 
transform. Apply your result to the transform you found in part (a). 



2.22 Compute the inverse transform, /(A:), for each of the following transforms: 



b) F(z) = 



z2-i.25z + 0.25' 



\z\ > 1; 




li 
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c) F(z) = 



z 



W>i; 



2.23 Use the z-transform to solve the difference equation 

y{k) - 3y{k - 1) + 2y{k - 2) = 2u(k - 1) - 2u(ifc - 2), 

{*, k>0, 
0, k<0, 
y(k) =0, ik < 0. 



CHAPTEF 

Sampl 



m 



3.1 INTR 

The use of di 
a continuous, 
pling. Sample 
tion, velocity, 
to calculate 1 
where discret 
in other parts 
continuous da 
of a purely co: 
is that of the 1 
is done via dij 
in the contini 
there as well, ' 
to discrete am 

quire some car 
be well reward' 
In this cha 
describe both ; 
also describe tl 
ing to recover 
analysis, we shi 



CHAPTERS 

Sampled-Data Systems 



3.1 INTRODUCTION 

The use of digital logic or digital computers to calculate a control action for 
a continuous, dynamic system introduces the fundamental operation of sam- 
pling. Samples are taken from the continuous, physical signals such as posi- 
tion, velocity, or temperature, and these samples are used in the computer 
to calculate the controls to be applied. Such digital controls are; hybrids, 
where discrete signals appear in some places and continuous signals occur 
in other parts. Such systems are called sampled-data systems because some 
continuous data are sampled before being used. In many ways, the analysis 
of a purely continuous system or of a purely discrete system is simpler than 
is that of the hybrid case. However, in digital control much of the processing 
is done via digital logic on discrete signals, but the origin of the signals is 
in the continuous world, and the destination of our computed outputs is 
there as well. Thus the role of sampling and the conversion from continuous 
to discrete and back from discrete to continuous are very important to im- 

P^t^RI^^Snl^ffl 

quire some careful treatment via the Fourier transform, but the effort will 
be well rewarded with the understandhig it will bring to later systems. 

In this chapter, we introduce the analysis of the sampling process and 
<tecribe both a time domain and a frequency domain representation. We 
jbo describe the companion process, that of sample extrapolation or hold- 
tag to recover a continuous time signal from its samples. As part of this 
ttialysis, we show that, because a sanapl^d-data system is made to be timp 
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varying by the introduction of sampUng, it is not possible to describe su^ 
eLtly by a transfer function. However, aRer sampling and holding} 
aCn^nuoS sfgnal is recovered, and we can approximate the response of, 
Imnle and hold to a sinusoid by fitting another sinusoid of the same ^ 
Tn y to the complete response. We wiU show how to compute this b.t-fi( 
sSdal respons; analyticaUy and experiment^y and thi^ have a good 
rrnnximation for a transfer function. For those familiar with the idea, our 
^"^ui^lent to the use of the "describing function" tl.t is used to^ 
app^oLate a transfer function for simple nonlinear systems. This concept: 

will be studied in Chapter 11. j ^ j 

Once the operations of sampUng and holding are miderstood, we wJl 
show that we can always represent the relationship between the sample o' 
thHnput and the samples of the output of a Unear constant ^tem by 
d^^l transfer function. Thus if we are wUUng to focus on the sample 
oSrthe entire power of Unear, constant system theory is available o us h 
Zt Z analysis of discrete linear systems is in many ways snnpler than tfe 
2ysl oStinuous linear systems in the way that subtraction is simpler; 

^^iation. We wiU see this when in the next c^^ we use 
tZ^f^m to analyze difference equations; however inthe digi ai con^ro| 
of^o^S^ dynamical systems we must understand the transitions from^ 
conUnuo^ to dbcrete and back again from discrete to continuous signals 
^^art. This is what we do in this chapter for the most elementary^ 
operations: sampling and holding. 

3.2 ANALYSIS OF THE SAMPLE AND HOLD 

.et samples of a physical signal such as a position or a velocity into digi^X 
fo'm weTjS^Uy have a sensor that produces a voltage Proport'-^^^'^ 
Z^kT^abl/and an analog-tcdigital converter, commmdy caUed^ 

conv^er or ADC, that transforms the voltage ii^o a 
xLpl^sical conversion always takes a fimte ^^^^ Jf^^^^^^ 
this Ume is significant with respect to the sample tune of ^h* co°<f ^^^^l 
^ ^^^t rSie rate of change of the signal being --P^^f^ 

is sketched in Fig. 3.1, where the switch, 5 is ^nelectromc device dri 
simple logic from a clock. Its operation is described betow. ^ 

With the switch, 5. in position I. the ampMer output ^out(^ fo^ « 
tracks the input voltage v^t) through the transfer function l/{RCs+lh 
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H = Hold; S in position 2. 
T = Track; 5 in position I. 



Figure 3.1 Analog-to-digital converter with sample and hold. 



circuit bandwidth of the SHC, 1/RC, is selected to be high compared to the 
input signal bandwidth. Typical values are iJ = 1000 ohms, C = 30 x 10 ^ 
farads for a bandwidth of / = l/27ri?C = 5.3 MHz. During this "trackmg 
time," the ADC is tm-ned oflF and ignores Vout- When a sample is to be taken 
at < = fcr the switch 5 is set to position 2 and the capacitor C holds the 
output of the operational amplifier frozen from that time at VontikT). The 
ADC is now signaled to begin conversion and has the constant input from 
the SHC to work on, so the resulting digital niunber is a true representation 
of the input voltage at the sample time. When the conversion is completed, 
the digfel mmiibgr is presented to tAie digital computer at which time the 
calculations based on this sample value can begin. The SHC switch is now 
moved to position 1, and the circuit is again tracking, waiting for the next 
conunand to freeze a sample. For example, the conversion time of the Bmr- 
Brown ADC803 is 1.5 microseconds for 12 bits of accuracy. The SHC needs 
only to hold the voltage for this short time m order for the conversion to 
be completed before it is started tracking again. The value taken is held 
inside the computer for the entire sampling period of the system, so the 
combination of the electronic SHC plus the ADC operate as a sample-and- 
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&^ voltag.^Mt!:^^^'^^^^^^^ - « t- representation 
«»f digital number is presented to fJi%^. f ^ conversion is completed, 
«fculations based on fhTs^]! * w**'^'t' '^P^t^^ which time the 
«»ved to position 1, l^^Ttl^^^"" The SHC switch is now 

command to freeze a ^^lZtZx:Vr^'^' ^^'""^ *^^'«-t 
Bwwn ADC803 is 1.5 microsec^nrr loV <^<>"ve«ion time of the Burr- 
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.^oompleted before it is st^t^t ' m order for the conversion to 
the computer for the^?^ ^^ain. The value taken is held 
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no 



k = -oo 



Figure 3.2 The sampler, 
hdd for the sampling period, T seconds, which may be many milliseconds 

For the purpose of the analysis we separate the sample and hold into 
two mathematical operations: a sampling operation represented as impulse 
modulation and a hold operation represented as a linear filter. The svmbd 
or schematic of the ideal sampler is shown m Fig. 3.2; its role is to rive a 
mathematical representation of the process of taking periodic samples from 

the sampled signals m the analysis of continuous signals using the Laplace 
transform.^ The technique is to use impulse modulation as the LthemTti^ 
representation of sampling. Thus, from Fig. 3.2, we picture the outputTf 
the sampler as a string of impulses. 



00 



'•*(0= E r(t)S{t-kT). 



(3.1) 



The impulse can be visualized as the limit of a pulse of unit area that has 
growmg amphtude and shrinkmg duration. The essential property of the 
impulse IS the.sifting property that oi tne 

f^jmt-a)dt = f{a) (3.2) 
the itTep"*"^ ^ continuous at o. The integral of the impulse is 

fj{r)dr^l(t), (3.3) 
and the Laplace transform of the unit impulse is one, because 



JWe assmne tlmt the reader has some femiliarity with Fourier and Laplace 
form analysis. For a general reference, see Bracewell (1978). 



(3.4) 
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?ourier and Laplace tra 
78). 



• E these properties we can see that r*(t), defined in (3.1), depends only 
the discrete sample values r(Jbr). The Laplace transform of r*(t) can be 
imputed as follows: 

£{r*(i)} = r r*iT)e-^^dT. 

J—OO 

If we substitute (3.1) for r*{t), we get 

= r T, riT)6{T-kT)e-'^dT, 



too ^ 

and now, exchanging integration and summation and using (3.2), we have 



(3.5) 



the notation R*{s) is used to symbolize the (Laplace) transform of r*{t), 
the sampled or impulse-modulated r(*).^ 

Having a model of the sampling operation as impulse modulation, to 
complete the description of the physical sample-and-hold we need to model 



H% will be necessary, from time to time, to consijder sampling a signal that is not 
continuous. The only case we will consider will be equivalent to applying a step 
function, to a sampler. For the purposes of this book we will define the unit 
atep to be continuous from the right and assume that the impulse, S{t), picks up 
the fiill value of unity. By this convention and (3.1) we compute 



1^(0 = E*(*-*^)' 



Jb=0 



and, using (3.2), we obtain 



(a) 



(b) 



I The reader should be warned that the Fourier integral converges to the average value 
?f a function at a discontinuity and not the value approached from the right as we 
assume. Because our use of the transform theory is elementary and the convenience 
I of equation (b) above is substantial, we have selected the continuous-from-the- 
,Hght convention. In case of doubt, the discontinuous term should be separated and 
tecated by special analysis, perhaps in the time domain. 



■f 
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3.5 



the hold operation, which will take the impulses that are produced by the! 
mathematical sampler and produce the piecewise constant output of the 
vice. Typical signals are sketched on Fig. 3.3. Once the samples are taken, as - 
represented by r*{t) in (3.1), the hold is defined as the means whereby these ^ 
impulses are extrapolated to the piecewise constant signal r^, defined as 



The requ: 



rh{t) = r(jkr) kT<t<kT + T. 



(3.6) i 



Because rh is composed of zero-order polynomials passed through the sam- 
ples of r{kT)j this hold operation is called the zero-order hold or ZOH and 
has the transfer function ZOH{s) . We can compute ZOH{s) by determining 
its impulse response. The hold filter will receive one unit-size impulse if the 
input signal is zero at every sample time except t — 0 and is equal to one 
there. In that case, r*(t) = 6{t) and r/|(t), which is now the impulse response 
of ZOH, is a pulse of height 1 and duration T seconds. The mathematical 
representation of the impulse response is, using the unit step function, 

p(0 = i(<)-i(<-r). 





Figure 3.3 The sample and hold, showing typical signals, (a) Input signal r; (b) 
sampled signal r*; (c) output signal r^; (d) sample and hold. 
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3.3 SPECTRUM OP A SAMPLED SIGNAL AND ALIASING 
The required transfer function is the Laplace transform of p(0 as 
ZOH(s) = C{p(t)} 

^fJ'iUt)-lit-T)]e-^dt 

Zl^^^y?^ 'TwI 'r^Z A/D converter with sample and hold can be 
Fig. 3.3 isLTfx^ted to rTpr^^^^^^^ -Snal r-(.) in 

.o^.ofthesamp.andJd^:^^^^ 



(3.7) 



J r; (b) 



AMn?r';''''''^ ^ SAMPLED SIGNAL 
AND ALIASING 

*=-<» n=-oo 

where the Fourier coefficipnf« r- • ,^ , 

^ coemcients, C„. are given by the integral over one period 

TLrn ^ - *r)e-M2T*/r)^< 
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bu. the sifting property ton. (M) this ea^ to integrate, with th. 

result 

Cn — ji* 

Thus we have derived the representation for the srun ot imputes as a Fourier 

series: 

(3.8) 



fc=-oo 



as the radian sampling frequency and now sub- 
(tsTJo r3.1) 'Is r « take the Up.«e transform o< the 
output of the mathematical sampler, 

we integrate the sum, term by term 

n=-oo ^ 

and if we combine the exponentials, 

n=-oo 

f ^f^^#^ Tn communication or radio en^neering 
where BW is the °™ f T^.g^t^^Z^ train corresponds to an 

ro^deiii^i'^J^sr::::?^'— 

in Fig. 3.4. 
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Figure 3.4 (a) Sketch of a spectrum amplitude and (b) the components of the 
spectnun after sampling, showing aliasing. 



An important feature of sampling, shown in Fig. 3.4, is illustrated at the 
frequency marked uji. Two curves are drawn comprising two of the elements 
of the sum at u;i. One of these, the larger amplitude in the figure located at 
the frequency ui, is the value of R{juJi)> The smaller component at ui comes 
from the spectrum centered at 2ir/T and is /2(ju;o), where uq is such that 
LJo-oJi- 2ir/T. This frequency, ldq, which shows up at o^i after sampling, 
is called in the trade an "alias" of a;i; the phenomenon is called aliasing. 

The phenomenon of aliasing has a clear meaning in time. Two contin- 
uous sinus6idi3 of different frequencies appear at the same frequency when 
sampled. We cannot, therefore, distinguish between them based on their 
samples alone. Fig. 3.5 shows a plot of a sinusoid at § Hz and of a sinusoid 
at g Hz. If we sample these waves at 1 Hz, as indicated by the dots, then 
we get the same sample values from both signals and would continue to get 
the same sample values for all time. Note that the sampling frequency is 1, 
and, if /i = i, then 
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Figure 3.6 (a) Sketch of a spectrum amplitude and (b) the components of the 
spectrum after sampling, showing removal of aliasing with an antialiasing filter. 



Band-limited signals that have no components above the Nyquist frequency 
are represented unambiguously by their samples. 

A corollary to the aliasing issue is the sampling theorem. We have seen 
that if R{jij) has componeftts above the Nyquist frequency u;s/2 or tt/T, then 
overlap and aliasing will occur. Conversely, we noticed that if R{juj) is zero 
for |u;| > tt/T, then sampling at intervals of T sec will produce no aliasing 
and the original spectrum can be recovered exactly from ii*, the spectrum 
of the samples. Once the spectrum is recovered, by inverse transform we can 
calculate the original sipiisfl itself. This is the sampling theorem: One can 
recover a signal from its samples, if the sampling frequency (a;, = 2n/T) 
is at least twice the highest frequency (tt/T) in the signal Notice that the 
sampUng theorem requires that R{juj) is exactly zero for all frequencies 
above n/T, 

A phenomenon somewhat related to aUasing is that of hidden oscilla- 
tions. If signal frequencies only up to tt/T* can be sampled without confusion, 
there is the possibility that a signal could contain some frequencies that the 
samples do not show at all Such signals, when they occur in a digital control, 
are called "hidden oscillations," an example of which is shown in Fig. 5.29. 
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3.4 DATA EXTRAPOLATION AND 
IMPOSTORS 

The sampUng theorem states that under the right conditions, it is possible 
to recover a signal from its samples; we now consider a formula for doii^ 

^h.J!''"l ^P^*^"™ «f ^(j^) i« contained in 

the low.frequency part o{R*{ju,). Therefore, to recover Vc) we need on v 
process R^iJu) through a low-pass filter and multiply by T to ^e^n R M 
a matter of fact, if RiJu,) has zero energy for frequencies'in the ba'n<^ fbot 
n/T {such an ^ is said to be band-Umited), then an ideal low-pass filter with 
gam r for - < a; < ./T and zero elsewhere would recover R{ju)lZ 

fr£T Tr^- l"^^T *^ '^'^^^ '°^-P^ characteristic I 

i^yju)). ihen we have the result 

R(Jw) = L{ju,)R*(Ju,). (3 

ffrL'^i!^/^*^ ^ the inverse transform of R{ju>), and because by (3.10) 
^^rl t^nsforms, its inverse transform r{t) must be 

he convolution of the time functions £(f) and r*it). The form of the filte 
unpulse response can be computed by using the definition of LUu,) from 
which the inverse transform gives^ ' 



1 r^/T 

J-kIT 



n/T 

T 



1 . TTf 

. vi 
= smc — . 
T 



(3.11) 



Using (3.1) for r*(<) and (3.11) for e{t), we find that their convolution is 
»-(0=rr(r) f; Sir - kT)sinc^^^^dT. 



^"Sinc'' is the name given to the function defined by smc{e) = sm{0)/9. 



0.8 



0.6 



0.4 



0.2 



-0.2 



-0.4 »— 



0^ 



Figured. 

Using the siftii 



Equation (3.12 
shows explicitl3 
r{t) from its sai 
time gaps betw 
Aplot of the in 
from the formu 

Because C*) is 
follows that thi 
starts at f = - 
t = 0! In man; 
needed until we 
be overcome by 
to the filter an« 
systems, a large 
approximations 



r 



3.4 DATA EXTRAPOLATION AND IMPOSTORS 113 




Figure 3.7 Plot of the impulse response of the ideal low-pass filter. 



Using the sifting property of the impulse, we have 

r(<)= f: r(A:r)sinc ^^^ I*^^^ 

ife=-oo 



(3.12) 



Equation (3.12) is a constructive statement of the sampling theorem: It 
shows explicitly how to construct the (by assumption) band-limited function 
r{t) from its samples. The sine fuiactions are the interpolators that fill in the 
time gaps between samples with a wave that has no frequencies above tt/T. 
A plot of the impulse response of this "ideal" hold filter is drawn in Fig. 3.7 
from the formula <tf (3.11). 

There is one senous (&Svir^^ the extrapolating signal given by (3.11). 
Because i{t) is the impulse response of the ideal low-pass filter L{ju)), it 
follows that this filter is noncausal because i{t) is nonzero for i < 0. l{t) 
starts at t = — oo when the impulse that triggers it does not occur until 
f = 0! In many communications problems the interpolated signal is not 
needed imtil well after the samples are acquired, and the noncausality can 
be overcome by adding a phase lag, e'^'^^y to L(jcj), which adds a delay 
to the filter and to the signals processed through it. In feedback control 
systems, a large delay is usually disastrous for stability, so we avoid such 
approximations to this function and use something else, like the polynomial 
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I'u^'^'rfi'''*^? zero-order hold already mentioned in connection with 
the ADC IS the most elementary. ^ 

hi Section 3.2 we introduced the zero-order hold as a model for thp 
storage register in an A/D converter that mamtains a constant signal valu. 
between samples. We showed m (3.7) that it has the transfer function 

1 — e"-'"^ 

ZOHUu;) = . (3 

We can discover the frequency properties of ZOH by expressing (3 13) 
m niagnitude and phase form. To do this, we factor out e--'-r/2 and multiDlv 
and divide by 2j to write the transfer function in the form 

The term in brackets is recognized as the sme, so this can be written 

wr/2 

and, using the definition of the smc function, 

ZOH{ju) = e~j'^'^f^Tsaic{uTl2). (3.14) 

Thus the effect of the zero-order hold is to introduce a phase shift of a;T/2 
corresponding to a time delay of T/2 seconds, and to multiply the gain by 
a function with the magnitude of 8inc(a;r/2). A plot of the magnitude L 
shown m Fig. 3.8, which iUustrates the fact that although the Lo^order 
hold IS a low-pass filter, it has a cut-off frequency weU beyond the Nyquist 
frequency The magnitude function is 



\ZOH(Ju,)\ = T\smc'^ 



(3.15) 



which slowly gets smaUer as u gets larger untU it is zero for the first tune 
at w = = 2?r/r. The phase is 

,lZOH{ju)) = ^^^^ (3.16) 
plus the 180» shifts where the sine function changes sign. 
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We can now give a complete analysis of the sample-and-hold circuit 
of Fig. 3.3 for a sinusoidal input r{t) in both the time and the frequency 
domains. We consider first the time domain, which is simpler, being just 
an exercise in construction. For purposes of illustration, we will use r{t) = 
38in(50^ + 7r/6) as plotted in Fig. 3,9. 

If we sample r{t) at the instants kT where the sampling frequency is 
w, = 2n/T = 2007r and T = 0.01, then the plot of the resulting rh{kT) is 
as shown in Fig. 3.9. Notice that although the input is a single sinusoid, the 
output is clearly not sinusoidal Thus it is not possible to describe this system 
ly a transfer function, because the fimdamental property of linear, constant 
^^fstsems is that a sinusoid iiiput produces sxi output^that. is a sinusoid of the 
same frequency, and the complex ratio of the amplitudes and phases is the 
transfer function. The sample-and-hold system is linear but tune varying. 
In the frequency domain, it is clear that the output r^it) contains more 
wan one frequency, and a complete analysis requires that we compute the 
amplitudes and phases of them all. However, in the application to control 
Bystems, the output of the hold will typically be applied to a dynamical 
^tem that is of low-pass character; thus the most important component in 
is the fundamental harmonic, at Uo = 50 rad/sec in this case. Also, in 
«te important field of digital filtering, one is usually using the digital filter to 
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Figure 3.9 Plot of 3sin(50t +7r /3) and the output of a sample-and-hold with samp 
period T = 0.01. 



replace an analog filter, and it is the fundamental that represents the sig 
component in the output. The other harmonics that appear in the outp 
are impostors, posing as signal when they are really unwanted consequen 
of the sample-and-hold process. In any event, we can proceed to analy 
rh{t) for all its harmonics and select out the fundamental component wh 
it makes sense, either by analysis or, in the implementation, by a low- 
anti-impostor filter. 

First, we need the spectrum of r(<). Because a simisoid can be d 
composed into two exponentials, we consider the Fourier transform 
v{t) = e^**'o*+"'*. For this we have 



J— oo 



(3.1 



This integral does not converge in any obvious way, but we can approach 
from the back door, as it were. Consider again the impulse, 6(1). The dir 
transform of this object is easy, considering the sifting property, as folio 
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If we exchange , with a, the integil nLdT ' «l>»n»tial 
Bluation 3.17 is of this form 

/oo 

Atth '^''^''("-.-o). ,3,jj 

to^i^X''^';'.^"i7ra' " to the delta 

°'^-ege„„and„„,tipBed,^?/^Tro?ri^„.:';i^^^ 
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Figure 3.10 Plot of the spectra of (a)/l, {h)R*, and (c)71a 



3.10(b) normalized by T. FinaUy, to find the spectrum of R^, we need only 
multiply the spectrum of /?• by the transfer function of the ^"^7^ 

ZOHUiv) = Te-i'^/hmc{o/r/2). 

Thus the spectriim of i?^ Is also a sum of an infinite number of terms, but 
now with intensities modified by the sine function and phases shiftJS^ the 
delay function u,T/2. These intensities are plotted in Fig. 3.10(c). NatS 

rnfadle^'t^T'"^ ".f " ^^'^^^'^-^ 

e^t^to Fig Sy*™ ^^'^^ Plotted 

nnlJlr ''^tu ^'^PP'^oximation to using only one sinusoid, we need 
Th^ ± o'^' ^^e first or fundamental harmonic from the componeite of i^^ 
This component has phase shift <f, and ampUtude A smc(a;T/2). In the time 
domain, the corresponding sinusoid is given by 

viit) = >l[smc(t.;r/2))sin[a;„(f - f)]. 

A plot of this approxunation for the signal from Fig 3.9 is given in Fig 3 11 

^lliit oflheT "^"^ sampled'and-held out^^t^f sho^ 

me nature of the approximation. 
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Figure 3.11 Plot of the output of the sample and hold and the first harmonic 
approximation. 



In control design, we can frequently achieve a satisfactory design for a 
sampled-data system by apprbximatmg the sample and hold with a con- 
tmuous transfer function corresponding to the delay of r/2. The controUer 
design is then done in the continuous domain but is implemented by com- 
puting a discrete equivalent. More discussion of this technique, sometimes 
caUed emulation, will be given in Chapter 5, where some examples iUustrate 
the results. 

Fiipsft-Order Hold. A more complex cu-cuit that might be preferred to 
the zero-order hold is the first-order hold, which, as suggested by its name, 
extrapolates data between sampling periods by a first-order polynomial,' 
a straight line. We can, as with the zero-order hold, compute the transfer 
function of the filter that, acting with impulse sampUng, produces the action 
of a first-order hold. A sketch of the response to a unit pulse will be helpful 
and is shown m Fig. 3.12. The first line, rising from 1 to 2 over the period 
from^ = 0 to i = r, is an extrapolation of the Ime between the points 
atj ^ -r (where the sample was, by assumption, equal to zero) and at 
t - 0, where the sample was unity. Likewise, the line going in the negative 
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direction from 0 at « = T to -1 at < = 2T is the extrapolation of th .• ' 
between the points at « = 0 and * = T. The Laplace transform of Il'^'l 
IS *'"is h[t)^ 



Obviously 
will be design* 
of improved p 



FOH(s) = 



(3.19)1 



The magnitude and phase of FOH{ji^) were plotted in Fig. 3.8 to I 
comparison with the characteristics of ZOH. Note that for low fre^.f^'^^l 
below n/2T) the first-order hold has significantly less phase lag t^^^Tl 
the ^ero-order hoW^ However, the FOH has much more ampUtude dTs w 
than the ZOH does. No clear guideHnes seem to exist inScati^ tt^^^^^^^ 
IS preferred to the other from the standpoint of an ideal contfol^sv '^^^ 
response; however, the increased hardware complexity of the FQIT ,1 
mentation almost always dictates that the ZOH be used. In soS^^f 
(e.g., hydraulic systems), the steps from the ZOH have been founr? ^ 
de^imental to the control actuator, and the solution has been to si^^ 
add a low-pass filter to the ZOH output with a time constant on thj o^dS 
of the sample period. FoUowing our earlier definition, this would be 3 
an antwmposter filter; some prefer the term "smoothing filter." S^ J 
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Figure 3.12 Impulse response of first-order-hold filter. 
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Figure 3.13 A cascade of samplers and filters. 



E*{s)G(s), then by definition we have 



but E*{s) is 



so that 



1 °° 

n=-oo 



k=-oo 



Now in (3.23) we can let ib = ^ - n to get 



<=-oo 



(3.22) 



(3.23) 



(3.24) 



oth« ''"jrds, because E* is already periodic, shifting it an integral number 
of perjQ^ leaves It un(^^, Su^iituiing, (.^.24) into (3.22) ySws 



U%s)^E'{s)l,Y^G{s-jru.,) 

— OO 

= £^*(«)(?*(5). 



(3.25) 



Note especially what is not true. If U(s) = E(<,\r(A tu^r. tt*( \ ^ 
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Figure 3.14 Block diagram of digital control as a sampled-data system. 



The final result we require is that, given a sampled-signal transform 
such as U*{s), we can find the corresponding z-transform simply by letting 
e*^ = 2 or 



(3.26) 



There is an important time-domain reflection of (3.26). The inverse Laplace 
transform of U*{$) is the sequence of impukes with intensities given by 
u{kT)\ the inverse z- transform of U{z) is the sequence of values u{kT), 
Conceptually, sequences of values and the corresponding z-transforms are 
easy to think about as being processed by a computer program, whereas the 
model of sampUng as a sequence of impulses is what allows us to analyze 
a discrete system embedded in a continuous world. Of course, the impulse 
modulator must always be followed by a low-pass circuit (hold circuit) in 
the physical world. Note that (3.26) can also be used in the other direction 
to obtain U*{s)y the Laplace transform of the train of impulses, from a given 
U{z). 

These rules of analysis can be illustrated by example. Consider the block 
diagram given in Fig. 3.14 taken from Fig, 1.1. In Fig. 3.14 we have modeled 
the A/D converter plus computer program plus D/A converter as an impulse 
modulator [which takes the samples from e(t)], a computer program that 
iprocesses t^hese samples, and a zero^prder hold that constructs the piecewise, 
constant output of the D/A corrverter from the impulses 6f trt*. In ttoefmtual 
computer we assume that the samples of e{t) are manipulated by a difference 
equation whose input-output effect is described by the z-transform D{z), 
These operations are represented in Fig. 3.14 as if they were performed 
on impulses, and hence the transfer function is D*{s) according to (3.26). 
Finally, the manipulated impulses, m*(t), are applied to the zero-order hold 
from which the piecewise-constant-control signal u{t) comes. In reality, of 
course, the computer operates on the sample values of e(t) and the piecewise- 
constant output is generated via a storage register and a D/A converter. The 
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impulses provide us with a convenient, consistent, and eflfective model of i 
processes to which Laplace-transform methods can be applied. < 
FVom the results given thus far, we caui write relations among Laplace 
trainsforms as 



E{s) = R-Y, 
M*{s) = E*D\ 



(3.27) 
(3.28) 

(3.29) 

(3.30) 

The usual idea is to relate the discrete output, V*, to the discrete input, fl*. 
Suppose we sample each of these equations by using the results of (3.5) to 
"star" each transform. The equations are^ 



Y = GU. 



E* = R*- Y*, 
M* = E*D*, 
U* = M*, 
Y* = [GU]\ 



(3.31) 
(3.32) 
(3.33) 
(3.34) 



Now (3.34) indicates that we need U, not U*, to compute Y*, so we must 
back up to, substitute (3.29) into (3.34): 



(3.35) 



Takuig out the periodic parts, which are those in which s appears only as 
e*^ [which include M*{s)], we have 



Substituting from (3.28) for M* gives 



(3^36) 



(3.37) 



^In sampling (3.29) we obtain (3.33) by use of the convention given in footnote 2, 
which follows (3.5) for impulse modulation of discontinuous functions. FVom the 
time-domain operation of the zero-order hold, it is clear that the samples of u and 
m are the same, and then from this (3.33) follows. 
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: substituting (3.31) for E* yields 

y* = (l-e-'^'')D*iG/8r[R*-Y*]. 



^Ifwe call 



then we can solve (3.38) for Y*, obtaining 



1 + 



(3.38) 



(3.39) 



"^^'^^ih^ (3.40) 
that the computer program corresponds to a discrete integrator 

^m = uikT-T) + Koe{kn (3.41) 

"rl'c^ttirS^^^^^ ''^'f - that the 
period T sotLte-"rTr\r^''^ Suppose we select the sampUng 
men in (TLl Pr.r~,^' "^"^^ *° ^^^^'"P"*^ components of 
function of (3.41), which m teri^ of z 

£?(^) -1-^-1 = jn- 

Using (3.26). we get the Laplace-transform form 



(3.42) 
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For the plant and zero-order-hold we require^ 

(i-.-'-)(c(.)/.r=(i-.---)(;5^^)" 

^ '\8 s + aJ 

Using (3.5), we have 

(1 - c-^')(G(.)/.r = (1 - e-^') {yI^s - i.e-ire-Ta ) • 
Because we assumed (for simplicity) that e"*^ = 5, this reduces to 



(l-e-^^)(G(.)/.)*==^i^^ 
_ 1/2 



1/2 



(3.43) 



(3.44) 



Combining (3.43) and (3.42), then, in this case, we obtain 

^ 2 (e*^-l)(e*^--l/2)* 

Equation (3.44) can now be used in (3.39) to find the closed-loop 
transfer function from which the dynamic and static responses can 
be studied, as a function of Kq, the program gain. We note also that 
beginning with (3.27), we can readily calculate that 



Y{s) = 



D* (1 - e-^') 



1 + iI* 



(3.45) 



Equation (3.45) shows how to compute the response of this system 
in between sampling instants. For a given r(t), the starred terms 
in (3.45) and the (1 - e^^*)-term correspond to a traun of impulses 
whose individual values can be computed by expanding in powers 
of e"^''. These impulses are applied to G{s)/s, which is the step 
response of the plant. Thus, between sampUng instants, we will see 
segments of the plant step response. 
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^Notice the similarity with (2.39) and Example 2.8. 
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Figure 3.15 A simple system that does not have a transfer function. 



With the exception of the odd-looking forward transfer function, (3.39) 
like the familiar feedback formula: forward-over-one-plus-feedback. Un- 
Ijrtunately, the sequence of equations by which (3.39) was computed was 
a bit haphazard, and such an effort might not always succeed. Another 
q t^wip le will further illustrate the problem. 



Example 3.2: Consider the block diagram of Fig. 3.15, which has 
only one sampling operation. This situation can arise if the error 
sensor has significant dynamics that precede the sampling action of 
the A/D converter. In this case H(a) represents the sensor dynamics. 
Again, we write the equations (all symbols are Laplace transforms): 



and again we sample 



E^R-Y, 
U = HE, 



E* = R*- Y*, 
U* = (HE)*, 



(3.46) 
(3.47) 
(3.48) 



(3.49) 
(3.50) 
(3.51) 



How do^we solve? In (3.50) we need E, not E\ So we must go back 

ir = {HiR-Y))* 
= (HR)* - (HY)*. 
Using (3.48) for Y, we have 

U*={HRy-{HU*G)\ 



Ill 
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Taking out the periodic U* in the second term on the right gives 
U*={HRy -U\HGy. 

Solving, we get 



i + {HGy 



(3.52) 



FVom (3.45), we can solve for Y* 



1 + {HG)*^ • 



(3.53) 



Equation (3.53) displays a curious fact. The transform of the input is 
bound up with His) and cannot be divided out to give a transfer functioni 
This system displays an important fact that all our facile manipulations 
of samples, and so on, might cause us to neglect: A sampled-data system 
IS ttme varying. The response depends on the time relative to the samvlina 
instant which the signal is applied. Only when the input samples alone are 
required to generate the output samples can we obtain a discrete transfer 
function. The time variation occurs on the takmg of samples. In general 
as m Fig. 3.15, the entire input signal r(t) is involved in the system re-' 
sponse, and the transfer-function concept fails. Even in the absence of a 
transfer function, however, the techniques developed here permit study of 
stability and response to specific inputs such as step, ramp, or sinusoidal 
Signals. 

We need to know the general rules of block-diagram analysis. In solving 
7 ^T** ourselves working with U, the signal that was sampled. 

f«wk#ss|m«»rMm^i^^ 

always "come free" after the equation sampling process and give a set of 
starred variables for which we can solve. 
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Example 3.3: Consider as a final example the block diagram 
drawn in Fig. 3.16. We select ^ and M as independent variables 
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Figure 3.16 A final example for transfer-function analysis of sampled-data 
terns. 



sjrs- 



and write 



M(s) = E*HGi. 



(3.54) 
(3.55) 



^dffrlm'(321> "se the "if periodic, then out" 



(3.56) 
(3.57) 



(3 Wiy ^^"^ equations by substituting for M* in (3.56) from 



E* = W-E*{HGiYGl 
R* 

l + {HGiyG^' 
lb obtain y we use the equation 

Y = E*H 

l + (jr<?i)*G5' 



(3.58) 



and 



l + iHGi)*Gi' 



(3.59) 



(3.60) 



this we have a transfer function. Why? Because only the 
samples of the external input are used to cause the output. To obtain 
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the z-transform of the samples of the output, we would let = z 
in (3.60). From (3.59) we can solve for the continuous output, which 
consists of impulses applied to ll{s) in this case. 



3.6 SUMMARY 

In this chapter we have considered the analysis of mixed systems that are 
partly discrete and partly continuous, taking the contmuous point of view. 
We used impulse modulation to represent the sampling process, and we 
derived the transfer functions of filters that would represent zero-order and 
first-order hold action. We showed that the transform of a sampled signal is 
periodic and that sampling introduces aliasing, which may be interpreted in 
both the frequency and the time domains. From the condition of no aliasmg 
we derived the sampling theorem. 

Fmally, we presented the block-diagram analysis of sampled-data sys- 
tems, showing that proper techniques, including the treatment of the sam- 
pler inputs as luiknowns, would lead to solution for the output transforms. 
However, we also found that not every sampled-data system has a transfer 
function. 



PROBLEMS AND EXERCISES 

3.1 Derive (3.45). 

3.2 Sketch a signal that shows hidden oscillations. 

3.3 Consider the circuit of Fig. 3.17. By plotting the response to a signal that is 
zero for all sample instants except i = 0 and that is 1.0 at « = 0, show that this 
circuit implements a first-order hold. 

3.4 Sketch the step response y(t) of the system shown in Fig. 3.18 for fc = |, 1, 
and 2. 




Figure 3.17 Block diagram of a first-order hold. 
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Figure 3.18 A sampled-data system. 



Figure 3.19 A general hold circuit. 

3.5 Sketcli the r^ponse of a second-order hold circuit to a step unit. What mieht 
be the major disadvantage of this data extrapolator? See Fig. 3 19 ^ 

8.6 Find the transform of the output Y(s) and its samples YUa) for the block 
^ shown m Fig. 3.20. Indicate whether a transfer'functi^ ^t^t 

3.7 Assume the foUowmg transfer functions are preceded by a zero-ord*.r hnH 
compute the resultmg discrete transfer functions zero-order hold and 

a) <?i(«) = l/a* 

b) G2(a) = €-i«'/(a + i) 



H 



(a) 



Hi. 



H 



(b) 



j^«-20 Block diagrams Of sampled data 



systems, (a) Single loop; (b) multiple 
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Y 






Gis) 






y(t) 

Figure 3.21 Block diagrams showing the modified z-transform. 

c) G3(s) = l/s{s + l) 

d) G7(s) = e-»»Va(« + 1) 

e) Gs{s) = l/{s^-l) 

3*8 One technique for examining the response of a sampled-data system between 
sampling instants is to shift the response a fraction of a sample period to the left 
and sample the result. The effect is as shown in the block diagram of Fig. 3.21 and 
is described by the equation 

m) = ir (5)Z{G(5)e"»^'}. 
The function Z{0{s)e"'^'} is called the modified z-transform of G{s). Let 



a) 



b) 



Compute y{t) by constructing the samples y{kT) from Y*{s) and ob- 
serving that with this plant, y{t) is an exponential decay with unit time 
constant over the ititersample interval. Sketch the response for five sample 
intervals. 

Let m = I and compute the samples correspondmg to Y*{s;m) (or 
y(z;m)]. Plot these on the same sketch as the samples of part (a) and 
verify that the midway points have been found. 

3.9 Show how to construct a signal of "hidden oscillations^" even one that grows 
in an unstable fashion but whose sample values are zero. Where in the ^-plane are 
the poles of the transforms of your signal(s)? 
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